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INTRODUCTION 


Computer  simulations  for  complete  saflng  and  arming  (S&A)  mechanisms  were 
developed.  Figures  1  and  2  show  the  essentials  of  two  different  configurations 
of  such  mechanisms,  consisting  of  spin  driven  rotors,  three  pass  involute  step-up 
gear  trains,  and  pin  pallet  runaway  escapements.  Figures  3  and  4  illustrate  two 
configurations  of  similar  mechanisms  which  contain  two  pass  involute  step-up  gear 
trains . 

The  present  work  draws  on  ideas  developed  previously  by  the  authors.  A 
dynamic  model  of  a  pin  pallet  runaway  escapement  driven  by  a  constant  torque  was 
developed  first  (ref  1).  The  ideas  of  pivot  and  tooth-to-tooth  friction  and  the 
resultant  etficiencies  of  multi-mesh  gear  trains  were  then  formulated  (ref  2), 
and  subsequently,  the  computational  technique  of  describing  mesh  initialization 
at  arbitrary  points  was  derived  (ref  3). 

This  report  reflects  many  of  the  ideas  originated  in  references  I  through 
3,*  in  describing  combined  mechanism  trains  which  are  driven  by  the  variable 
torques  of  centrifugal  rotors.  Again,  the  three  motion  regimes  of  the  runaway 
escapement,  i.e.,  coupled  motion,  free  motion,  and  impact,  are  considered.  In 
addition,  the  ability  to  include  the  effect  of  a  pallet  which  has  an  arbitrarily 
located  center  of  mass  has  been  introduced.  Further,  all  non-impact  contact 
forces  can  now  be  determined  for  strength  considerations. 

The  complete  and  detailed  derivations  of  both  mathematical  models  and  their 
configuration  variations  are  given.  The  associated  computer  programs  are  de¬ 
scribed  in  detail,  with  instructions  for  running  them.  These  programs  were  ex¬ 
tensively  tested,  both  for  the  configuration  no.  1  data  (fig.  1)  of  the  three 
pass  M125A1  booster  S&A  as  well  as  for  the  configuration  no.  2  data  (fig.  4)  of 
the  two  pass  M577  SSD  mechanism.  The  M125A1  booster  simulation  armed  in  4b  turns 
for  average  friction  conditions  at  various  spin  rates.  Under  similar  conditions, 
the  simulation  of  the  M577  SSD  required  31.7  turns  for  its  arming  cycle. 


DESCRIPTION  OF  COMPUTER  PROGRAMS  SANDA3  AND  SANDA2 


The  essential  steps  of  programs  SANDA3  and  SANDA2 ,  which  are  listed  in  ap¬ 
pendixes  C  and  D,  respectively,  follow.  This  discussion  is  based  on  the  computer 
program  description  in  reference  1,  which  deals  with  a  runaway  escapement  driven 
by  a  constant  torque.  For  the  sake  of  clarity  the  description  in  reference  1  is 
repeated  and,  whenever  necessary,  changes  in  the  programs  due  to  the  presence  of 
the  gear  trains  and  the  rotors  are  identified. 
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These  models  can  be  adapted  to  spring  driven  mechanisms  with  only  minor  changes. 


f|N,0N  3  .  N.  X  FUZE  BODY 


Figure  1.  Rotor  driven  S§A  device  with  three  pass  step-up  gear  train 
configuration  no.  1 


S5A  device  with  two  pass  step-up  gear  train. 


Both  programs  may  be  used  for  each  of  the  two  possible  mechanisms  configu¬ 
rations.  To  differentiate  between  these  configurations,  a  special  input  para¬ 
meter  CONFIG  must  be  provided. 

The  flow  chart  of  the  pin  pallet  runaway  escapement  part  of  the  programs,  is 
shown  in  figure  5  which  was  taken  from  reference  1.  The  choice  of  variable  des¬ 
ignations  was  made  so  that  they  differ  as  little  as  possible  from  the  nomencla¬ 
ture  used  in  the  various  derivations  given  in  the  appendixes  of  the  present  re¬ 
port  as  well  as  in  appendixes  A  through  F  and  H  of  reference  1 . 

Both  main  programs  start  with  the  reading  in  and  writing  of  the  configura¬ 
tion  type  as  well  as  of  all  relevant  physical  data.  This  is  followed  by  the 
computation  of  gear  ratios  (app  A  and  B),  fuze  body  angles  of  the  applicable 
configuration,  with  the  help  of  signum  function  s^  (app  A  and  B  and  ref  2), 
centrifugal  forces  (ref  2),  earliest  and  latest  possible  values  of  the  gear  an¬ 
gles  by  way  of  subroutine  ALFA  (ref  2),  and  the  initialization  of  the  gear  angles 
(ref  3).  In  each  case,  the  simulation  begins  with  entrance  coupled  motion  at  a 
starting  angle  PHIL),  which  represents  the  angle  <j>  of  the  escape  wheel  which  is 
associated  with  the  approximate  center  of  the  escape  wheel  tooth.  This  angle 
corresponds  to  the  cumulative  escape  wheel  angle  PHITOT  of  0°. 


Program  SANDA3 

Coupled  Motion  (Location  1) 

To  solve  the  differential  equation  of  coupled  motion  (eq  A-18S,  app  A), 
the  main  program  calls  on  an  available  fourth  order  Runge-Kutta  routine.  The 
main  purpose  of  the  subroutine  FCT  is  to  present  the  second  order  differential 
equation  in  terms  of  two  first  order  equations  to  RKGS.  PHI(l)  and  PHI(2)  repre¬ 
sent  the  angle  <j>  and  the  angular  velocity  $  ,  respectively.  The  computation  of 
all  parameters  of  the  differential  equation  takes  place  in  subroutine  FCT  as  well 
as  in  subroutines  IN3  and  KINEM. 

_  •  • 

The  subroutine  KINEM  computes  current  values  of  g,  g,  i(j  and  includ¬ 

ing  the  moment  arms  AJ,  B^,  Cj,  and  Dj  (which  are  primed  in  the  present 
report).  The  associated  expressions  for  the  above  values  are  cited  in  reference 

1. 


Subroutine  IN3  computes  various  gear  mesh  parameters  and  instantaneous 
mesh  contact  angles,  as  well  as  the  signum  functions  Sj  to  s5.  In  addition,  the 
parameters  Aj  to  A53,  listed  in  appendix  A,  are  obtained.  (Note  that  the  program 
uses  the  symbols  AA1 ,  etc.,  and  that  the  parameters  A  54  to  A  5^  are  computed  in 
subroutine  FCT.  This  division  is  dictated  by  the  limited  number  of  arguments 
permitted  in  any  one  subroutine.)  Finally,  the  gear  indexing  operation  (ref  2) 
is  performed  with  the  help  of  the  angle  4. 


^  RKGS  Routine,  IBM  System/360  Scientific  Subroutine  Package,  (360A-CM-0X3) 
Version  III. 
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Figure  5.  Flow  chart  for  simulation  of  pin  pallet  runaway  escapement 


The  manner  in  which  the  instantaneous  rotor  angle,  $  +  N,,^,  of  the 

coupled  motion  differential  equation  (A-185)  must  be  expressed  in  subroutine  FCT 
will  now  be  discussed.  Recall  that  $  is  the  initial  rotor  angle,  ^  is  the 
total  angle  of  rotation  of  the  escape  wheel,  and  stands  for  the  gear  ratio 
between  the  rotor  and  the  escape  wheel.  Since  the  angle  <(>,  and  consequently  the 
Runge-Kutta  variable  PHI(l),  varies  between  140°  and  155°  during  entrance  and 
between  185°  to  200°  during  exit  coupled  motion,  the  total  escape  wheel  angle 
4>  can  only  be  obtained  by  continuously  adding  the  increments  due  to  each  cycle 
of  Runge-Kutta  computations.  Thus 


where 


V)T 


Aip 


\0T  +  A* 


total  escape  wheel  angle  up  to  a  given 
Runge-Kutta  cycle.  This  is  represented 
by  PHITOT  in  the  program. 

increment  of  escape  wheel  during  a  given  Runge- 
Kutta  cycle. 


This  increment,  A<p,  is  calculated  as  the  difference  between  the  latest  value  of 
FH 1(1)  and  its  previous  value,  which  has  been  stored  as  PHIPR.  In  this  manner, 
equation  1  becomes 


<hT 


PHITOT  +  PHI(l)  -  PHIPR 


(2) 


Subroutine  FCT  also  decides  on  the  values  of  I  and  1^  as  required  by 
equations  A-88  and  A-89  as  well  as  equations  A- 155  and  A-156,  respectively.  The 
associated  conditional  statements  assign  the  larger  values  for  these  combined 
moments  of  inertia  whenever  the  product  of  the  angular  velocity  and  the  angular 
acceleration  is  positive,  i.e.,  both  quantities  have  the  same  sign  (app  F  and 
paragraph  preceding  eq  A-148). 

The  associated  subroutine  OUTP  is  responsible  for  printing  out  the 
results  <fr,  $,  and  $  together  with  the  current  values  of  time,  g,  g,  <j»,  and 
PHITOT.  Further,  all  coupled  motion  contact  forces  are  calculated  according  to 
equations  A-192  to  A-196,  and  the  maximum  values  of  these  forces  during  one  rotor 
cycle  are  determined.  Finally,  this  subroutine  decides,  according  to  the  follow¬ 
ing  criteria,  whether  coupled  motion  is  to  be  continued.  These  are: 

1.  The  distance  g  <  0  (ref  1,  eq  A-12)  and  because  of  the  nature 
of  the  coordinate  system  g  is  always  negative  while  the  pallet  pin  makes  contact 
with  the  escape  wheel  tooth,  and 

2.  The  contact  force  P  >0  (eq  A-196) 

n 
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The  introduction  of  the  contact  force  between  the  escape  wheel  and  the  pallet 
aLlows  the  use  of  simpler  controls  than  those  given  in  reference  1. 

Thus,  when  control  is  returned  to  the  main  program,  it  is  either  be¬ 
cause  the  pallet  pin  has  left  the  end  of  the  tooth  and  there  is  no  further  possi¬ 
bility  of  coupled  motion,  or  the  pin  has  disengaged  from  the  inside  of  the 

tooth.  In  either  case,  free  motion  results  and  control  is  eventually  shifted  to 
the  subroutines  FCTF  and  OUTPF  through  RKGS  (location  5).  This  is  done  directly 
if  g  <  0.  If  g  >  0,  the  main  program  must  decide  whether  the  preceding  computa¬ 
tions  have  been  made  for  entrance  or  for  exit  action  and  whether  the  next  contact 
will  occur  in  enrance  or  in  exit.  In  the  sample  mechanism,  g  =  0  when  <p  is 
approximately  155°  for  entrance  action  and  approximately  200°  for  exit  action. 
If  g  >  0  and  $  <  160°,  all  possibility  for  entrance  contact  is  ended  and  (J>  must 
be  incremented  by  the  tooth  angle  S  (ref  1,  figs.  2  and  5),  while  ip  must  be  in¬ 
cremented  by  the  angle  (2  if  -  X)  in  order  to  get  the  initial  angles  for  exit 
motion.  For  g  >  0  and  <}>  >  160°,  entrance  contact  is  expected  and  <(>  must  be 

decremented  by  the  angle  2  S  (ref  1,  fig  4,  where  the  new  entrance  tooth  no.  3 
comes  into  action).  At  the  same  time,  the  pallet  angle  ip  must  be  decremented  by 
( —2 tt  +  A)  for  the  initial  angle  of  entrance  motion.  This  shift  from  one  side  of 
the  pallet  to  the  other  must  also  be  observed  in  connection  with  the  pallet  cen¬ 
ter  of  mass  angle  (figs.  1  and  2).  During  entrance  action, 


*c  -  PS  ICC  (3) 

in  the  program,  while  during  exit  action 

♦c  -  A  +  PSICC  (4) 

These  indexing  operations  due  to  pallet  position  occur  both  in  coupled  and  in 
free  motion.  They  have  no  effect  on  the  continuous  computation  of  the  cumulative 
escape  wheel  angle  PHITOT. 


Free  Motion  (Location  5) 

The  differential  equations  of  free  motion,  as  given  by  equation  A-197 
for  the  pallet  and  equation  A-203  for  the  combined  escape  wheel,  gear  train  and 
rotor,  are  again  solved  by  the  Runge-Kutta  routine.  In  order  to  obtain  the  mag¬ 
nitudes  of  the  variables  4>  and  ^  and  their  derivatives  at  identical  times,  the 
two  independent  second-order  differential  equations  are  transformed  into  four 
simultaneous  first-order  equations.  (While  only  the  two  first-order  equations 
associated  with  each  of  the  two  variables  are  actually  coupled,  the  routine 
treats  all  four  as  if  they  were  coupled  and  thus  produces  solutions  for  identical 
time  increments.)  These  four  expressions,  which  are  presented  in  subroutine  FCTF 
are  given  in  the  following  form. 


DX( 1 )  -  X(2)  :  (-  |) 


(5) 
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DX(3) 


X(4)  : 


(6) 


(-  *) 

DX(2)  =  - —  [-A65(X(2))2  +  Ab6X(2)  +  Ab/ 

+  A68stn(  4>1C  +  N4l(<(>T  +  X(l)  -  PHIPR) )  ] 

:  (=  +)  (7) 

DX(4)  =  ~  [-A21(X(4))2  -  A61X(4)  -  Ab2 

"60 

+  A63sin(Yp-i!)  -  4»c>]  :  (=  ijO  (8) 


Subroutine  FCTF  also  computes  the  parameters  Agg  to  Agg 
IN3  for  the  computation  of  all  gear  related  parameters. 
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and  calls  on  subroutine 


The  associated  subroutine  OUTPF  computes  the  free  motion  contact  forces 
according  to  equations  A-209  through  A-211,  and  finds  their  maxima,  after  deter¬ 
mining  the  parameters  A^q  to  Agg.  The  computation  of  these  latter  variables  also 
requires  a  call  on  subroutine  IN3.  In  addition  to  the  above,  a  continuous  count 
of  PHITOT  is  provided  in  OUTPF.  This  angle,  time,  <j>,  <£,  ij,  and  the  contact 
forces  are  printed.  This  subroutine  also  makes  the  decision  of  whether  or  not  to 
remain  in  free  motion.  The  sensing  variables,  f  and  g  "  “  GP  (ref  1,  eq  D-7  and 
D-Ll),  are  used  for  this  purpose. 

If  f  >  0  and  g'  <  0,  free  motion  is  continued  without  indexing. 
If  f  >  0  and  g'  >  0,  free  motion  is  also  continued,  but  since  contact  is  no 
longer  possible  for  the  component  pair  for  which  the  previous  computations  were 
made,  indexing  must  take  place.  This  is  accomplished  in  the  same  manner  as  in 
coupled  motion. 


If  f  <0,  control  is  returned  unconditionally  to  the  main  program.  If 
it  finds  that  g'  >0,  indexing  takes  place  and  control  is  given  back  to  the  sub¬ 
routine  FCTF.  If  g'  <  0,  contact  is  about  to  take  place  or  has  just  occurred. 
This  program  must  decide  whether  this  contact  represents  just  a  close  approach 
which  is  followed  by  further  free  motion,  an  impending  impact,  or  the  beginning 
of  coupled  motion.  To  this  end,  the  quantities  V  and  Vg  (ref  1,  eq  F-22  and  F- 
23),  are  computed  for  the  entrance  and  exit  free  motion  tests. ^  The  first  three 


3 

Whenever  A  ■  Ip^  *  0,  the  simulation  stops  because  of  division  by  zero. 

Should  this  occur,  FCTF  prints  "IPR  EQUALS  ZERO  -  SIMULATION  TERMINATED." 

^  Under  the  present  circumstances,  if  $  <  160°,  only  entrance  contact  can  follow; 
if  $  >  160°,  exit  contact  will  occur. 
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cases  of  the  entrance  free  motion  test  of  the  main  program  are  illustrated  in 
figure  6.  With  both  angular  velocities  (  =  PHI(2)  =  X(2)  and  p  -  DPSI  -  X(4)) 

positive,  the  following  three  possibilities  exist: 

1 •  If  |  Vp|  >  j  Vs|  ,  the  contacting  surfaces  will  separate  again, 
free  motion  will  result,  and  control  must  be  transferred  to  subroutine  RJCGS  (lo¬ 
cation  6),  which  calls  free  motion  subroutine  FCTF. 

2.  If 

let  in  coupled  motion, 
tion  1  initiates  RKGS),  which  in  this  case,  calls  coupled  motion  subroutine  FCT. 

3  •  If  |  V  |  <  J  V  (  ,  impact  will  occur,  and  control  must  be  given 

to  subroutine  IMPACT  (location  15). 


|  V  |  -  |VS|  ,  the  escape  wheel  will  start  driving  the  pal- 
and  control  must  be  transferred  to  subroutine  RKGS  (loca- 


.  »  • 

(  and  are  positive. 


Distance  f  is  enlarged.) 


Figure  6.  Entrance  free  motion  test 


Impact  (Location  15) 

.  ,  The  subroutine  IMPACT  uses  the  current  values  of  the  .angular  velocities 

and  ^  and  computes  the  post  impact  angular  velocities  and  4>  ,  applying 
equations  F— 20  and  F— 21,  ref  1.  (The  moment  of  inertia  of  the  escape  wheel  is 
now  expressed  according  to  equation  A-212,  appendix  A  of  the  present  report, 
which  refers  the  rotor,  as  well  as  the  gear  train  inertias,  to  the  escape  wheel 
shaft.  As  shown  in  reference  1,  the  tangential  impact  has  been  neglected  and, 
therefore,  E2  *  and  F2  ■  Aj .) 

After  control  is  returned  to  the  main  program,  it  is  decided  whether 
free  or  coupled  motion  follows  the  impact.  This  is  accomplished  by  considering 
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the  post-impact  contact  point  velocities,  Vp  and  Vg,  in  the  impact  tests,  which 
are  similar  to  the  free  motion  tests. 

If  the  contact  velocities  are  vectorially  equal  to  each  ether,  or  if 
the  absolute  value  of  the  difference  of  their  absolute  magnitudes  is  less  than 
2.0  in. /sec  (5.08  cm/sec),  control  is  transferred  to  coupled  motion.  If  these 
velocities  indicate  a  subsequent  separation,  which  is  more  usual,  computation  is 
transferred  to  free  motion. 

Reversal  of  Gear  Train  Motion  Due  to  Impact 

If  the  impact  torque  on  the  escape  wheel  is  sufficiently  large,  the 
motion  of  the  gear  train  may  be  temporarily  reversed  (i.e.,  the  escape  wheel 
angular  velocity  |  becomes  negative).  This  causes  the  friction  forces  between 
the  gear  teeth  and  at  the  various  gear  pivots  to  be  reversed  in  direction.  (The 
normal  forces  between  the  gear  teeth  remain  unaffected,  and  the  normal  bearing 
forces  are  solved  for  in  the  usual  manner.)  This  change  in  the  direction  of  the 
friction  forces  is  expressed  both  for  coupled  and  free  motion  by  letting  the 
coefficient  of  friction  u  of  all  gear  train  components  become  negative  (app  E) . 
This  is  accomplished  in  subroutines  IN3  and  IN2  by  the  following  use  of  the  sig- 
num  function  |/  |  : 


MU  =  ABS(MU)*  |/  ||| 


(9) 


[The  coefficient  of  friction  associated  with  the  escapement  interface  and  the 
pallet  pivot  is  called  pj  (read  into  the  program  as  MU1).]  Any  motion  reversal 
at  these  locations  is  accounted  for  by  the  signum  functions  s^  and  ,  respec¬ 
tively. 

Termination  of  Computations 

Computations  are  terminated  whenever  the  geared  motion  of  the  rotor 
ends.  This  corresponds  to  $  =  PHICUTD.  (See  Program  SANDA3,  General  Parameters, 
below.)  The  duration  of  the  subsequent  unretarded  motion  of  the  rotor  is  assumed 
to  be  negligible. 


Program  SANDA2 

Program  SANDA2  is  very  similar  to  program  SANDA3,  and  where  possible,  the 
following  description  will  refer  to  program  SANDA3. 

Coupled  Motion  (Location  L) 

To  solve  the  differential  equation  of  coupled  motion  (app  B,  eq  B-113), 
the  main  program  calls  on  the  fourth  order  Runge-Kutta  routine.  Subroutine  FCT 
presents  the  second  order  differential  equation  in  terms  of  two  first  order  equa¬ 
tions  to  RKGS,  which  computes  PHI(l)  and  PHI(2)  as  in  SANDA3 .  The  computation  of 
all  necessary  parameters  takes  place  in  FCT  and  KINEM.  In  addition,  subroutine 
IN2  is  used.  This  subroutine,  which  differs  from  IN3,  computes  various  gear  mesh 
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parameters  and  contact  angles  of  the  two  pass  train.  It  further  computes  the 
signum  functions  S|,  s^,  and  85  as  well  as  the  parameters  A^  to  A^,  listed 

in  appendix  B.  The  parameters  A45  to  Ajq  are  computed  in  subroutine  FCT.  Final¬ 
ly,  the  gear  indexing  operation  is  performed  with  the  help  of  the  angle 

The  instantaneous  rotor  angle  now  becomes  and  is  treated 

in  FCT  in  the  same  manner  as  shown  by  equations  1  ana  2.  Subroutine  FCT  also 
decides  on  the  values  of  IpR  and  1^  as  required  by  equations  A-88  and  A-89  as 
well  as  equations  B-87  and  B-88.  The  associated  conditional  statements  assign 
the  larger  values  for  these  combined  moments  of  inertia  whenever  the  product  of 
the  angular  velocity  and  the  angular  acceleration  is  positive,  i.e.,  both  quanti¬ 
ties  have  the  same  sign.  (See  discussion  following  eq  2.) 

The  associated  subroutine  OUT?  is  responsible  for  printing  out  the 
results  <(>,  I,  and  $,  together  with  the  current  values  of  time,  g,  g,  ip,  iji,  and 
PHITOT.  Further,  all  coupled  motion  contact  forces  and  their  maxima  are  deter¬ 
mined  according  to  equations  B-120  to  B-123.  The  coupled  motion  control  criteria 
given  in  this  subroutine  are  identical  to  those  described  in  connection  with 
subroutine  OUTP  of  SANDA3.  This  is  also  true  of  the  indexing  operations  con¬ 
nected  with  the  angles  <}»,  <|f,  and  i|>c. 

Free  Motion  (Location  5) 

The  differential  equations  of  free  motion,  as  given  in  equation  B-124 
for  the  pallet  and  equation  B-130  for  the  combined  escape  wheel,  gear  train  and 
rotor,  are  also  solved  by  the  Runge-Kutta  routine.  The  associated  four  simulta¬ 
neous  first  order  differential  equations,  which  are  given  in  subroutine  FCTF,  are 
similar  to  those  used  in  SANDA3 .  (See  paragraph  before  eq  5.). 


DX(  1) 

=  X(2)  :  (=  |) 

(10) 

DX(3) 

11 

X 

/-N 

W 

/-v 

» 

-e* 

(11) 

DX(2) 

~  [-A56(X(2))2  -  A57X(2)  +  AbiJ 

«55 

-A598in(<J>lc  +  N31(<f>T  +  X(l)  -  P11IPR) )  ] 

:  (=  4>) 

(12) 

DX(4) 

“  ~  [-A21(X(4))2  -  A52X(4)  -  A 53 

A51 

+  A5l+sin(Yp  -  t|>  -  4>c)]  :  (*  ♦)  U3) 


The  subroutine  FCTF  also  computes  the  parameters  A^  to  A^^  and  calls  on  subrou¬ 
tine  IN2  for  the  computation  of  all  gear  related  parameters. 
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The  associated  subroutine  OUTPF  computes  the  parameters  A^j  to  Acg  in 
order  to  obtain  the  free  motion  contact  forces  according  to  equations  B-136  and 
B-137.  In  addition,  it  determines  the  maxima  of  these  contact  forces.  These 
operations  also  require  a  call  on  subroutine  IN2.  Again,  a  continuous  count  of 
PHITOT  is  provided  in  OUTPF,  and  time,  <j>,  $,  \J>,  ij>  and  the  contact  forces  are 
caused  to  be  printed  together  with  PHITOT.  Whether  or  not  to  remain  in  free 
motion  is  decided  by  the  same  free  motion  tests  outlined  in  SANDA3. 

Impact  (Location  15) 

The  treatment  of  the  impact  phase  of  the  motion  is  identical  to  that 
given  in  connection  with  SANDA3.  (The  moment  of  inertia  of  the  escape  wheel  is 
now  expressed  according  to  equation  B-138,  app  B,  which  refers  the  rotor  as  well 
as  gear  and  pinion  no.  2  to  the  escape  wheel  shaft.) 

Reversal  of  Gear  Train  Motion  Due  to  Impact 

Any  reversal  of  the  gear  train  due  to  impact  is  treated  in  the  same 
manner  as  that  described  in  SANDA3  (eq  9).  The  control  of  the  coefficient  of 
friction,  p,  now  takes  place  in  subroutine  IN2.  The  coefficient  of  friction 
associated  with  the  escapement  interface  and  the  pallet  pivot  is  again  called 
Pj  and  is  read  into  the  program  as  MU1. 

Termination  of  Computations 

Computations  are  terminated  whenever  the  geared  motion  of  the  rotor 
ends.  This  corresponds  to  <p  =  PHICUTD  (Program  SANDA2 ,  General  Parameters, 
below).  The  duration  of  the  subsequent  unretarded  motion  of  the  rotor  is  assumed 
to  be  negligible. 


COMPUTER  SIMULATIONS  OF  EXAMPLE  MECHANISMS 


Computer  simulations  of  modified  M125A1  and  M577  spin  driven  S&A  devices  are 
discussed  in  this  section.  Both  contain  appropriately  adapted  involute  gear 
trains  rather  than  the  actual  clock  gear  trains  used.-’  Program  SANDA3  has  been 
used  for  the  M125A1  mechanism,  which  has  configuration  type  no.  1.  Both  the 
program  as  well  as  the  relevant  output  are  given  in  appendix  C.  The  simulation 
of  the  M577  mechanism,  which  has  configuration  no.  2,  was  accomplished  with  the 
help  of  SANDA2.  This  program,  as  well  as  associated  output,  is  given  in  appendix 
D.  Both  simulations  were  run  for  30,000  RPM  in  order  to  obtain  maximum  contact 
forces. 


^  All  gears  and  pinions  have  the  original  number  of  teeth  and  diametral  pitches. 
The  individual  meshes,  with  unity  contact  ratio,  are  designed  to  operate 
at  standard  center  distances  and  have  unequal  addenda  (obtained  with  program 
INV0L1  of  ref  2). 
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The  data  requirements  of  both  programs,  an  explanation  of  the  various  output 
data,  and  the  manner  in  which  the  number  of  turns  to  arm  results  are  obtained  for  a 
given  spin  velocity  are  presented  below. 


Program  SANDA3 
Input  Data 

The  required  input  data  is  the  first  portion  of  the  output  of  Program 
SANDA3.  They  represent  the  mechanism  parameters  of  the  M125A1  booster  and  are 
listed  below  both  as  computer  variables  and  as  symbols. 


A 

B 

C 


Escapement  Parameters. 

a  =  0.2061  in.  (5.234  mm) 


distance  between  pivots  Op  and  0S  (fig.l) 


=  0.192  in.  (4.877  mm)  ■  escape  wheel  radius 
-  0.0788  in.  (2.001  mm) 


pallet  radius  from  pivot  to  center  of  pin 
(identical  for  entrance  and  exit) 


R 

r  =  0.015  in, 

.  (0.381  mm)  -  pallet  pin  radius  (identical  for 

entrance  and  exit) 

ALPHA  = 

0  =  51° 

*  escape  wheel  tooth  half  angle 

CONFIG  » 

1 

-  configuration  no.  1 

EREST  =» 

e  »  0 

■  coefficient  of  restitution  (high 
speed  motion  pictures  of  runaway 
escapement  suggest  totally 
inelastic  impacts) 

LAMBDA  - 

X  -  152.144“ 

«  angle  between  entrance  and  exit  pallet 
radii 

DELTA  - 

S  -  30° 

■  angle  between  individual  escape  wheel 
teeth 

For  details  of  the  above  nomenclature  see  reference  1. 

Mass  Properties  of  Components. 

Ml  “  mj 

-  0.1155  x  10"3 

lb-sec^/in.  (2.022  x  10“^  kg)  *  mass  of  rotor 

1'2  ■  m2 

-  0.686  x  10"5 

lb-sec^/in.  (1.201  x  10-3  kg)  -  mass  of  gear  and 

pinion  no.  2 

M3  ■  m3 

-  0.282  x  10"5 

lb-sec^/in.  (4.938  x  10“3  kg)  -  mass  of  gear 

and  pinion  no.  3 
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r 


r 


i 

i 


M4  =  m4  -  0.252  x  10"5  lb-sec2/in.  (4.412  x  10-4  kg) 


mass  of  escape 
wheel  and  pinion 
no.  4 


0.359  x  10~5  lb-sec2/in.  (6.286  x  10“4  kg)  =  mass  of  pallet 


MP  = 

11  - 

12  » 

13  = 

14  =* 


% 


X  =  0.1055  x  10"4  in. -lb-sec2  (0.1189  x  10"5  kg-m2)  =  moment  of 

1  <nprtla  o 


inertia  of 
rotor 


I0  *  0.126  x  10"6  in. -lb-sec2  (0.142  x  10*7  kg-m2)  =  moment  of 


inertia  of  gear 
and  pinion  no.  2 


I3  *>  0.308  x  10 


_7  7  ,rt-8  .  2)  =  moment  of  inertia 

7  in. -lb-sec^  (0.348  x  10  °  kg-m 


of  gear  and  pinion 
no.  3 


IP 


1,  -  0.268  x  10"7  in. -lb-sec2  (0.303  x  10"8  kg-m2)  -  moment  of 

^  of  escape  wheel 

and  pinion  no.  4 

t  =  0.245  x  10'7  in. -lb-sec2  (0.277  x  10“8  kg-m2)  =■  moment  of  inertia 

p  of  pallet 


General  Parameters. 


RCl 

* 

rCl 

m 

0.106 

in.  (2.692 

RCP 

as 

rCP 

m 

0 

RHOP 

** 

pp 

m 

0.016 

in.  (0.406 

RPM 

=* 

30,000 

PHI1CD 

s 

^lC 

- 

144° 

PSICCD 

m 

*c 

m 

0° 

PHID 

145° 

mpi)  -  distance  from  rotor  pivot  to  center 
of  mass 

*  pallet  eccentricity 
mm)  ■  pallet  pivot  radius 

*  spin  rate 

-  rotor  angle  in  starting  position 
(fig.  1) 

■  eccentricity  angle  of  pallet 
(fig.  1) 

-  escape  wheel  starting  angle 
of  initial  coupled  notion 
simulation  (for  choice  of 
this  angle,  see  ref  1) 
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PHICUTD 


3600 


MU 


MU1 


PSUBD1  - 

PSUBD2  - 

PSUBD3  - 

NG1 

NG2 

NG3 

NP2 

MP3 

NP4 

CAPRP1  - 

CAPRP2  - 


■  cumulative  escape  wheel 
angle,  obtained  from  product 
of  total  rotor  rotation 
and  gear  train  ratio.  The 
total  rotor  rotation  for  the 
M125A1  is  76°  while  the  gear 
ratio  is  47.25.  Thus,  PHICUTD 
-  76  x  47.25  -  3600°. 


u  =  0.15  =  coefficient  of  friction  of 

gear  train  (pivots  and  tooth-to- 
tooth  contacts,  includes  escape 
wheel  pivot) 

Ml  ”0.15  ”  coefficient  of  friction  on  pallet- 

escape  wheel  interface  and  pallet 
pivot 

Gear  Parameter. 


Pdl  ”  44 
Pd2  -  65 
Pd3  -  77 

NCi  -  42 


”  diametral  pitch  of  mesh  no.  1 
(rotor  and  pinion  no.  2) 

”  diametral  pitch  of  mesh  no.  2 
(gear  no.  2  and  pinion  no.  3) 

”  diametral  pitch  of  mesh  no.  3 
(gear  no.  3  and  escape  wheel 
pinion) 

■  number  of  teeth  of  rotor  (full 
gear  no.  1) 


nG2 

%3 

Np2 

Np3 

NP4 


-  27 
”  27 

-  8 

-  9 
”  9 


■  number  of  teeth  of  gear  no.  2 

■  number  of  teeth  of  gear  no.  3 

■  number  of  teeth  of  pinion  no. 

■  number  of  teeth  of  pinion  no. 

■  number  of  teeth  of  pinion  no. 

(escape  wheel  pinion) 


2 

3 

4 


Rpl  ”  0.47727  in.  (12.1226  mm)  ■  pitch  radius  of  gear  no.  1 
r  (rotor) 


Rp2  ■  0.20769  in.  (5.2753  mm)  ■  pitch  radius  of  gear  no.  2 
Rp3  *  0.17532  in.  (4.4531  mm)  *  pitch  radius  of  gear  no.  3 


CAPRP3 


RP2 

- 

rp2 

am 

0.09091  In.  (2.3091  mn) 

3 

pitch  radius  of  pinion  no.  2 

RP3 

m 

rp3 

= 

0.06923  In.  (1.7584  mm) 

- 

pitch  radius  of  pinion  no.  3 

RP4 

= 

rp4 

m 

0.05844  In.  (1.4844  mm) 

“ 

pitch  radius  of  pinion  no.  4 
(escape  wheel  pinion) 

THETA1 

=■ 

0i 

- 

20° 

= 

pressure  angle  of  mesh  no.  1 

THETA2 

m 

02 

=» 

20° 

= 

pressure  angle  of  mesh  no.  2 

THETA3 

- 

03 

3 

20° 

= 

pressure  angle  of  mesh  no.  3 

R1 

- 

«1 

3 

0.225  in.  (5.715  mm) 

3 

distance  of  rotor  pivot  from 
spin  axis 

R2 

=* 

6(  2 

■ 

0.436  in.  (11.074  mm) 

= 

distance  of  pivot  of  gear  and 
pinion  set  no.  2  from  spin  axis 

R3 

=x 

d?  ^ 

- 

0.504  in.  (12.602  mm) 

= 

distance  of  pivot  of  gear  and 
pinion  set  no.  3  from  spin  axis 

R4 

3M 

** 

0.520  in.  (13.208  mm) 

= 

distance  of  pivot  of  escape 
wheel  from  spin  axis 

R5 

m 

«5 

- 

0.433  in.  (10.998  mm) 

* 

distance  of  pivot  of  pallet 
from  spin  axis 

RH01 

m 

Pi 

- 

0.062  in.  (1.575  mm) 

= 

pivot  radius  of  rotor 

RH02 

m 

P2 

- 

0.025  in.  (0.635  mm) 

= 

pivot  radius  of  gear  and  pinion 
no.  2 

RH03 

m 

P3 

m 

0.018  in.  (0.457  mm) 

- 

pivot  radius  of  gear  and  pinion 
no.  3 

RH04 

m 

P4 

m 

0.016  in.  (0.406  mm) 

- 

pivot  radius  of  escape  wheel 

CAPRB1 

m 

Rbl 

m 

0.44849  in.  (11.3916  mm) 

3 

base  radius  of  gear  no.  1 

CAPRB2 

3 

Rb2 

- 

0.19517  in.  (4.9573  tmn) 

- 

base  radius  of  gear  no.  2 

CAPRB3 

am 

Rb3 

at 

0.16475  in.  (4.1847  an) 

B 

base  radius  of  gear  no.  3 

RB2 

am 

rb2 

m 

0.08543  in.  (2.1699  mm) 

m 

base  radius  of  pinion  no.  2 

RB3 

s 

rb3 

am 

0.06506  in.  (1.6525  an) 

m 

base  radius  of  pinion  no.  3 

RB4 

91 

rb4 

m 

0.05492  in.  (1.3950  mm) 

- 

base  radius  of  escape  wheel 
pinion 

CAPROl 

m 

Rol 

m 

0.48791  in.  (12.3929  mm) 

- 

outside  radius  of  gear  no.  1 
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CAPR02 

= 

Ro2 

= 

0.21579 

in. 

(5.4811  mm) 

a 

outside 

radius 

of 

gear  no.  2 

CAPR03 

= 

Ro3 

= 

0.18216 

in. 

(4.6269  ram) 

a 

outside 

radius 

of 

gear  no.  3 

R02 

= 

ro2 

= 

0.110  in 

1.  (2 

.7940  mm) 

a 

outside 

radius 

of 

pinion  no.  2 

R03 

= 

ro3 

= 

0.08089 

in. 

(2.0546  mm) 

a 

outside 

radius 

of 

pinion  no.  3 

R04 

= 

ro4 

3 

0.06828 

in  • 

(1.7343  mm) 

= 

outside 

pinion 

radius 

of 

escape  wheel 

J1 

3 

J1 

= 

0 

3 

initialization 

parameter  for  mesh 

no.  1.  (The  value  of  zero 
corresponds  to  earliest  possible 
contact  of  mesh,  ref  3.) 


J2  =  J2  =  0  =  initialization  parameter  for  mesh 

no .  2 

J3  =  J3  =  0  *  initialization  parameter  for  mesh 

no .  3 


Output  Data 

Fuze  Geometry.  The  angles  BETA1D  =  to  BETA4D  »  |34  and  GAMMA2D  =  Y2  to 
GAMMA5D  =  are  printed  for  checking  purposes. 

Coupled  motion.  For  each  time  increment  T  of  the  coupled  motion,  the 
following  variables  are  evaluated: 


PHI 

3 

4> 

3 

instantaneous  escape  wheel  angle  (deg) 

PH I DOT 

= 

• 

♦ 

S 

escape  wheel  angular  velocity  (rad/sec) 

G 

= 

g 

= 

pallet  -  escape  wheel  contact  position  (in.),  ref  1 

GDOT 

3 

• 

g 

a 

time  rate  of  change  of  g  used  to  determine  the 
direction  of  relative  sliding  (in. /sec),  ref  1 

PS  ID 

= 

<l> 

a 

pallet  angle  (deg) 

PSIDOT 

3S 

- 

pallet  angular  velocity  (rad/sec) 

PHITOT 

m 

\ 

a 

cumulative  escape  wheel  angle  (deg) 

F34 

a 

F34 

a 

normal  contact  force  of  gear  no.  3  on  pinion 
no.  4  (lb) 

F23 

m 

F23 

m 

normal  contact  force  of  gear  no.  2  on  pinion 
no.  3  (lb) 
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F12 

3 

F12  = 

normal  contact  force  of  gear  no.  1  on  pinion 
no.  2  (lb) 

PN 

= 

j-c 

II 

normal  contact  force  between  escape  wheel  and 
pallet  (lb),  computed  according  to  equation  A-195 

PNPSI 

II 

-3- 

C 

a* 

normal  contact  force  between  escape  wtieel  and 
pallet  (lb),  computed  according  to  equation  A-196 
(serves  for  checking) 

DPMI  2 

= 

♦  = 

escape  wheel  acceleration  ( rad/sec  ■0,  Runge- 
Kutta  output 

Free  Motion.  For  each  time  increment  T  of  the  free  motion,  the  following 
variables  are  evaluated: 


PHI 

<t>  = 

instantaneous  escape  wheel  angle  (deg) 

PH I DOT 

= 

£ 

escape  wheel  angular  velocity  (rad/sec) 

PSI 

= 

♦  = 

pallet  angle  (deg) 

PS  I  DOT 

= 

11 

pallet  angular  velocity  (rad/sec) 

PHITOT 

= 

^T 

cumulative  escape  wheel  angle  (deg) 

FF12 

= 

FF12  = 

normal  contact  force  of  gear  no.  1  on  pinion 
no.  2  for  free  motion  (lb) 

FF23 

= 

FF23  = 

normal  contact  force  of  gear  no.  2  on  pinion 
no.  3  for  free  motion  (lb) 

FF34 

= 

FF34  = 

normal  contact  force  of  gear  no.  3  on  escape 
wheel  pinion  for  free  motion  (lb) 

Impact .  Just  preceding  the  "IMPACT”  label,  the  program  prints  the  values 
VP  and  VS,  which  stand  for  the  pre-impact  velocity  components,  normal  to  the  escape 
wheel  tooth,  of  both  the  pallet  and  escape  wheel  contact  points  (ref  1).  Subsequent 
to  the  "IMPACT"  label,  the  following  variables  are  evaluated: 


PHI 

♦  = 

instantaneous  escape  wheel  angle  (deg),  same  as 
before  impact 

PHIDOT 

3 

• 

♦  = 

post-impact  escape  wheel  angular  velocity  (rad/sec) 

PSI 

- 

<t> 

pallet  angle  (deg),  same  as  before  impact 

PSIDOT 

= 

post-impact  pallet  angular  velocity  (rad/sec) 

PHITOT 

*S 

*r 

cumulative  escape  wheel  angle  (deg),  same  as 
before  impact 
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VP 


post-impact  normal  velocity  component  of  pallet 
(pin)  at  contact  point  (ref  1) 


VS  =  Vg  =  post-impact  normal  velocity  component  of  escape 

wheel  tooth  at  contact  point  (ref  1) 

In  the  present  program  VP  is  equal  to  VS  since  the  coefficient  of  restitution  is 
zero. 

Turns-to-Arm  and  Maximum  Contact  Forces.  The  number  of  turns-to-arm  is 
obtained  by  the  program  for  that  time  T^qq  which  corresponds  to  the  escape  wheel 
angle  PHIGUTD  =  3600°.  Therefore 

RPM 

Number  of  turns-to-arm  **  — —  X  T.j600 

Thus,  with  T-j^qq  ®  0.09068  second  (the  actual  time  corresponds  to  3600.04°)  and 
30,000  RPM 

Number  of  turns-to-arm  =  — -  X  0.09068 

oU 

=  45.34  turns 


The  maximum  non-impact  contact  forces  for  the  total  cycle,  both  for  coupled  and  for 
free  motion,  are  listed  at  the  end  of  the  output. 


Program  SANDA2 
Input  Data 

The  required  input  data  represent  the  first  portion  of  the  output  of 
program  SANDA2.  They  stand  for  the  following  mechanism  parameters  of  the  M577  SSD: 

Escapement  Parameters. 


A 

a 

=  0.2097  in.  (5.326  mm) 

=  distance  between  pivots  Op 
and  Og 

B 

= 

b 

=  0.193  in.  (4.902  mm) 

=  escape  wheel  radius 

C 

c 

=■  0.0788  in.  (2.001  mm) 

=  pallet  radius  from  pivot  to 
center  of  pin  (identical  for 
entrance  and  exit) 

R 

St 

r 

=  0.015  in.  (0.381  mm) 

=  pallet  pin  radius  (identical 
for  entrance  and  exit) 

ALPHA 

S3 

a 

**  52° 

=  escape  wheel  tooth  half 
angle 

r 


CONFIG 

2 

■  configuration  no.  2 

EREST 

e  =0 

«  coefficient  of  restitution 
(high  speed  motion  pictures 
of  runaway  escapement 
suggest  totally  inelastic 
impacts) 

LAMBDA 

= 

X  =  152.144° 

=  angle  between  entrance  and 
exit  pallet  radii 

DELTA 

= 

5  =  30° 

=  angle  between  individual 
escape  wheel  teeth 

Details  of 

the  above  nomenclature  are 

given  in  reference  1. 

Mass  Properties  of  Components. 

Ml  =  = 

0.3531 

X  10-4  lb-sec2/ln. 

(6.182  X  10“^  kg)  =  mass  of  rotor 

M2  =  = 

0.3226 

X  10-5  lb-sec2 /In. 

(5.648  X  10-4  kg)  =  mass  of  gear  and 

pinion  no.  2 

M3  =  m3  = 

0.2656 

X  10"5  lb-sec2 /In. 

(4.650  X  10-4  kg)  =  mass  of  escape 

wheel  and  pinion 
no.  3 

MP  =  mp  = 

0.2451 

X  10"5  lb-sec2 /in. 

(4.291  X  10-4  kg)  =  mass  of  pallet 

11  =  Ix  = 

0.3350 

X  10-5  in. -lb-sec2 

(0.3788  X  10”6  kg-m2>  =  momen,:  of  ln*rtia 

of  rotor 

12  =  I2  - 

0.3750 

X  10~7  in. -lb-sec2 

(0.4240  X  10-8  kg-tn2)  »  moment  of  inertia  1 

gear  and  pinion  no 

13  »  I3  = 

0.2490 

X  10-7  in. -lb-sec2 

(0.2815  X  10  8  kg-m2)  =  moment  of  inertia  1 

escape  wheel  and 
pinion  no.  3 

IP  =  Ip  =  °*6572  X  10-7  In. -lb-sec2 

General  Parameters. 

(0.7430  X  10-8  kg-m2)  =  moment  of  inertia  < 

pallet 

RC1 

rCl 

=  0.115  in.  (2.921 

nm)  ■  distance  from  rotor  pivot  to  center 
of  mass 

RCP 

rCP 

-  0 

■  pallet  eccentricity 
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RHOP  =  Py  =  0.0155  in.  (0.394  mm)  =*  pallet  pivot  radius 
RPM  =  30,000  ■  spin  rate 


PHI1CD  =  <j>  -  -122.23° 

PSICCD  =  t|>,  =0° 

c 

PHID  =  150° 

PHICUTD  -  1600° 

MU  =  u  »  0 . 1 5 

MU1  =  ui  0.15 


=  eccentricity  angle  of  rotor 
in  starting  position  (fig.  4) 

«  eccentricity  angle  of  pallet 
(fig.  4) 

=  escape  wheel  starting  angle  of 
initial  coupled  motion  simulation 
(for  choice  of  this  angle,  see 
pp  7  and  14,  ref  1) 

*  cumulative  escape  wheel  angle, 
obtained  from  product  of  total 
rotor  rotation  and  gear  train  ratio. 
The  total  rotor  rotation  for  the 
M577  is  48.3°  while  the  gear  ratio 
is  33.02.  Thus,  48.3  X  33.02  »  1600. 

=  coefficient  of  friction  of  gear 
train  (pivots  and  tooth-to-tooth 
contacts,  includes  escape  wheel 
pivot) 

=  coefficient  of  friction  on  pallet  - 
escape  wheel  interface  and  pallet 
pivot 


Gear  Parameters. 


PSUBD1 

ss 

pdl 

44 

s 

diametral  pitch  of  mesh  no.  1  (rotor 
and  pinion  no.  2) 

PSUBD2 

= 

Pd2 

s 

63.5 

= 

diametral  pitch  of  mesh  no.  2  (gear 
no.  2  and  escape  wheel  pinion) 

NG1 

- 

NG1 

= 

41 

= 

number  of  teeth  of  rotor  (full  gear 
no.  1) 

NG2 

s 

NG2 

m 

29 

- 

number  of  teeth  of  gear  no.  2 

NP2 

- 

Np2 

m 

6 

m 

number  of  teeth  of  pinion  no.  2 

NP3 

- 

Np3 

m 

6 

- 

number  of  teeth  of  pinion  no.  3 
(escape  wheel  pinion) 

CAPRP1 

- 

RP1 

m 

0.46591  in. 

(11.834  mm)  ■ 

pitch  radius  of  gear  no.  1  (rotor) 
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CAPRP2 

zm 

RP2 

-  0.22835 

in 

.  (5.800  mm) 

RP2 

=a 

rp2 

-  0.06818 

in 

.  (1.732  mm) 

RP3 

= 

rp3 

=  0.04724 

in 

.  (1.200  mm) 

THETA1 

SB 

01 

-  20° 

THETA2 

- 

02 

=  20“ 

R1 

- 

«1 

m  0.225  in. 

(5.715  mm) 

R2 

- 

H  2 

■  0.408  in. 

(10.200  mm) 

R3 

- 

<«3 

»  0.368  in. 

(9.100  mm) 

R4 

3 

=■  0.361  in. 

(9.025  mm) 

RH01 

3 

Pi 

-  0.0465 

in. 

(1.181  mm) 

RH02 

3 

P2 

-  0.0155 

in. 

(0.394  mm) 

RH03 

= 

P  3 

-  0.0155 

in. 

(0.394  mm) 

CAPRB1 

- 

Rbl 

-  0.43781 

in 

.  (11.120  nm) 

CAPRB2 

= 

Rb2 

**  0.21458 

in 

.  (5.450  mm) 

RB2 

= 

rb2 

-  0.06407 

in 

.  (1.627  mm) 

RB3 

3 

rb3 

-  0.04439 

in 

.  (1.128  mm) 

CAPROl 

= 

Rol 

-  0.47388 

in 

.  (12.037  mm) 

CAPR02 

- 

Ro2 

-  0.23387 

in 

.  (5.940  mm) 

R02 

S 

ro2 

=■  0.09373 

in 

.  (2.381  mm) 

R03 

3 

ro3 

-  0.07321 

in 

.  (1.860  mm) 

J1 

m 

Jl 

-  0 
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pitch  radius  of  gear  no.  2 

pitch  radius  of  pinion  no.  2 

pitch  radius  of  pinion  no.  3 
(escape  wheel  pinion) 

pressure  angle  of  mesh  no.  1 

pressure  angle  of  mesh  no.  2 

distance  of  rotor  pivot  from 
spin  axis 

distance  of  pivot  of  gear  and 
pinion  set  no.  2  from  spin  axis 

distance  of  pivot  of  escape  wheel 
from  spin  axis 

distance  of  pivot  of  pallet  from 
spin  axis 

pivot  radius  of  rotor 

pivot  radius  of  gear  and  pinion 
no.  2 

pivot  radius  of  escape  wheel 

base  radius  of  gear  no.  1 

base  radius  of  gear  no.  2 

base  radius  of  pinion  no.  2 

base  radius  of  escape  wheel  pinion 

outside  radius  of  gear  no.  1 

outside  radius  of  gear  no.  2 

outside  radius  of  pinion  no.  2 

outside  radius  of  escape  wheel 
pinion 

initialization  parameter  for  mesh 
no.  1.  (The  value  of  zero 
corresponds  to  earliest  possible 
contact  of  mesh  (ref  3).) 


J2  =*  J2  -  0 


»  initialization  parameter  for  mesh 
no.  2 


Output  Data 

Fuze  Geometry.  The  angles  BETA1D  ■  to  BETA3D  *  and  GAMMA2D  »  to 
GAMMA4D  *  are  printed  for  checking  purposes. 


Coupled  Motion.  For  each  time  increment,  T,  of  the  coupled  motion,  the 
following  variables  are  evaluated: 


PHI 

3 

ai 

instantaneous  escape  wheel  angle  (deg) 

PHIDOT 

- 

• 

♦ 

3 

escape  wheel  angular  velocity  (rad/sec) 

G 

= 

g 

3 

pallet  -  escape  wheel  contact  position  (in.)  (ref  1) 

GDOT 

- 

• 

g 

= 

time  rate  of  change  of  g  used  to  determine  the  direction 
of  relative  sliding  (in. /sec)  (ref  1) 

PSID 

= 

♦ 

= 

pallet  angle  (deg) 

PSIDOT 

3 

* 

3 

pallet  angular  velocity  (rad/sec) 

PHITOT 

= 

= 

cumulative  escape  wheel  angle  (deg) 

F23 

3 

F23 

- 

normal  contact  force  of  gear  no.  2  on  pinion  no.  3  (lb) 

F12 

SB 

F12 

= 

normal  contact  force  of  gear  no.  1  on  pinion  no.  2  (lb) 

PN 

= 

Pn 

8 

normal  contact  force  between  escape  wheel  and  pallet  (lb), 
computed  according  to  equation  B-83 

PNPSI 

- 

Pn^ 

= 

normal  contact  force  between  escape  wheel  and  pallet  (lb), 
computed  according  to  equation  B-84,  (serves  for  checking) 

DPHI2 

= 

♦ 

- 

escape  wheel  acceleration  (rad/sec  ),  Runge-Kutta  output 

Free  Motion.  For  each  time  increment,  T,  of  the  free  motion,  the  follow¬ 
ing  variables  are  evaluated: 


PHI  *  <J>  =  instantaneous  escape  wheel  angle  (deg) 

PHIDOT  *  |  =  escape  wheel  angular  velocity  (rad/sec) 

PSI  ■  <J»  *  pallet  angle  (deg) 
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PSIDOT 


♦ 


PHITOT 

$T 

FF12 

*  FF12 

FF23 

=  FF23 

=  pallet  angular  velocity  (rad/sec) 

=  cumulative  escape  wheel  angle  (deg) 

=  normal  contact  force  of  gear  no.  1  on  pinion  no.  2 
for  free  motion  (lb) 

=  normal  contact  force  of  gear  no.  2  on  escape  wheel 
pinion  for  free  motion  (lb) 


Impact.  Just  preceding  the  "IMPACT"  label,  the  program  prints  the  values 
VP  and  VS,  which  stand  for  the  pre-impact  velocity  components,  normal  to  the  escape 
wheel  tooth,  of  both  the  pallet  and  escape  wheel  contact  points  (ref  1).  Subsequent 
to  the  "IMPACT"  label,  the  following  variables  are  evaluated: 


PHI 

* 

ss 

instantaneous  escape  wheel  angle  (deg),  same  as  before 
impact 

PHID0T 

= 

• 

4> 

ae 

post-impact  escape  wheel  angular  velocity  (rad/ sec) 

PSI 

a 

♦ 

= 

pallet  angle  (deg),  same  as  before  impact 

PSIDOT 

a 

• 

♦ 

post-impact  pallet  angular  velocity  (rad/sec) 

PHITOT 

= 

4>t 

= 

cumulative  escape  wheel  angle  (deg),  same  as  before  impact 

VP 

= 

VP 

a 

post-impact  normal  velocity  component  of  pallet  (pin)  at 
contact  point  (ref  1) 

VS 

= 

Vs 

= 

post-impact  normal  velocity  component  of  escape  wheel  tooth 
at  contact  point  (ref  1) 

Note 

that  in 

the 

present 

program  VP  is  equal  to  VS  since  the  coefficient  of  restitu- 

tion 

is  zero 

• 

Turns-to-Arm  and  Maximum  Contact  Forces.  The  result  of  27 .93  turns  to  arm 
at  a  spin  velocity  of  30,000  RPM  is  obtained  by  the  program  in  a  manner  similar  to 
that  shown  in  Program  SANDA3  under  paragraph  entitled  Turns-to-Arm  and  Maximum 
Contact  Forces  for  the  MI25AI  booster  mechanism.  (Now  the  escape  wheel  angle 
PHICUTD  -  1600°.) 

The  maximum  non-impact  contact  forces  for  the  total  cycle  both  for  coupled 
and  free  motion,  are  listed  at  the  end  of  the  output. 
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CONCLUSIONS 


While  It  was  not  the  purpose  of  the  present  investigation  to  undertake  para¬ 
metric  studies  of  either  of  the  two  mechanisms  for  which  the  programs  were  written, 
both  programs  were  sufficiently  tested  to  confirm  that  such  studies  are  possible. 
They  may  include  variations  in  masses  and  moments  of  inertia  of  all  components, 
variations  of  gear,  escapement  and  fuze  geometries,  as  well  as  various  friction 
conditions. 

Test  runs  for  average  friction  conditions  and  standard  geometries  showed  that 
the  M577  SSD  simulation  required  27.9  turns  to  arm  regardless  of  the  assumed  spin 
rate.  Under  similar  conditions,  the  M125A1  booster  model  needed  45.3  turns  to 
arm.  Both  of  these  values  are  well  within  the  military  specifications  for  these 
devices . 

The  ability  to  determine  the  contact  forces  during  coupled  and  free  motion 
furnishes  the  foundation  for  future  work  concerning  the  strength  of  fuze  gear 
teeth.  Such  work  will  also  have  to  include  the  determination  of  the  appropriate 
dynamic  factors  due  to  the  impact  of  the  runaway  escapement. 

A  most  interesting  phenomenon,  which  may  also  serve  as  the  starting  point  of  a 
future  investigation,  occurred  following  the  first  exit  impact  during  the  initial 
motion  cycle  of  the  1-1577  SSD  simulation.  The  preceding  coupled  motion  was  initiated 
closer  to  the  root  of  the  escape  wheel  tooth  than  usual,  i.e.,  the  starting  angle  of 
the  escape  wheel  was  chosen  to  be  PHID  =  145°,  when  a  more  usual  value  during  stable 
operation  is  equal  to,  or  larger  than  150°.  With  this  entrance  starting  condition 
and  the  resulting  higher  than  usual  pallet  angular  velocity,  the  exit  impact  was 
sufficiently  large  to  reverse  the  motion  of  the  escape  wheel  and  the  gear  train. 
During  this  backward  motion  of  the  escape  wheel,  its  deceleration  (i.e.,  accelera¬ 
tion  in  the  forward  direction),  due  to  the  influence  of  the  driving  rotor,  became 
approximately  20  times  larger  than  usual  and  all  contact  forces  increased  substan¬ 
tially.  This  increase  in  forces  subsided  as  soon  as  the  forward  motion  of  the 
escape  wheel  was  restored.  This  phenomenon  did  not  occur  during  normal  operation, 
since  once  the  escapement  has  become  stable  coupled  motion  starts  much  nearer  to  the 
tip  of  the  escape  wheel  tooth.  The  resulting  smaller  free  motion  angular  velocities 
of  the  pallet  caused  less  severe  impacts  and,  as  a  result,  prevented  motion  rever¬ 
sals. 
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APPENDIX  A 

DYNAMICS  OF  ROTOR  DRIVEN  S&A  MECHANISM 
WITH  A  THREE-PASS  INVOLUTE  GEAR  TRAIN 
AND  A  PIN  PALLET  RUNAWAY  ESCAPEMENT 


\ 


J 


■! 


29 


DESCRIPTION  OF  SYSTEMS  AND  OUTLINE  OF  DERIVATIONS 


The  present  appendix  gives  derivations  for  a  complete  mathematical  model  of  an 
S&A  mechanism,  consisting  of  a  spin  driven  rotor,  a  three  pass  involute  step-up  gear 
train,  and  a  pin  pallet  runaway  escapement.  Figures  A-l  and  A-2  show  the  two  types 
of  configurations  to  which  this  model  can  be  adapted. 

This  work  draws  to  a  considerable  extent  on  previous  efforts  by  the  authors, 

i.e.,  the  dynamics  of  the  pin  pallet  runaway  escapement  (ref  A-l)  and  the  fuze  gear 
train  analysis  (ref  A-2).  As  in  reference  A-l  the  following  three  regimes  of  the 
mechanism  are  considered:^  * 

1 .  Coupled  Motion 

The  escape  wheel  is  in  constant  contact  with  one  of  the  pallet  pins  while 
it  is  driven  by  the  rotor  (gear  no.  1)  through  the  gear  and  pinion  sets  no.  2  and 
3.  The  coupled  motion  differential  equation  in  the  escape  wheel  variable  $  is 
obtained  by  combining  the  solutions  to  the  Newtonian  force  and  moment  expressions 
for  the  individual  mechanism  components. 

2.  Free  Motion 

The  pallet  and  the  escape  wheel  -  gear  train  -  rotor  systems  move  Inde¬ 
pendently  of  each  other  in  this  phase  of  motion.  Thus,  there  results  one  differen¬ 
tial  equation  for  the  pallet  in  the  variable  <l>  ,  and  another  one  for  the  combined 
system  in  the  escape  wheel  variable  if. 

3.  Impact 

The  formulation  of  the  impact  regime  is  taken  directly  from  reference  A-l, 
except  that  now  the  moment  of  inertia  of  the  escape  wheel  and  pinion  no.  4  also 
contains  the  referred  mass  properties  of  the  rotor  and  gear  and  pinion  sets  no.  2 
and  3.  This  impact  simulation  is  based  on  the  classical  angular  impulse  -  momentum 
model,  where  a  coefficient  of  restitution  is  used  to  account  for  the  energy 
losses.  Since  it  is  assumed  that  the  effect  of  the  impact  force  between  the  escape 
wheel  and  the  pallet  is  much  greater  than  the  effects  due  to  the  driving  torque  of 
the  rotor  and  the  various  retarding  torques  due  to  friction,  the  latter  are  disre¬ 
garded. 

The  influence  of  friction  forces  is  considered  both  in  the  coupled  and  the  free 
motion  regimes.  There  is  friction  at  the  escape  wheel  -  pallet  interface  during 
coupled  motion,  and  there  is  friction  between  the  gear  teeth  and  at  all  pivots 
during  both  of  these  regimes.  As  in  reference  2  the  individual  pivot  friction 


For  better  understanding,  consult  figures  in  reference  A-l. 
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gure  A-l.  Rotor  driven  SfjA  device  with  three  pass  step-up  gear  train  -  configuration  no. 


ee  pass  step- 


torques  are  obtained  by  the  algebraic  addition  of  the  two  friction  moments  due  to 

the  x  and  y  components  of  the  normal  bearing  forces,  rather  than  by  the  direct  use  ' 

of  the  resulting  normal  forces.  This  conservative  approach  to  friction  is  necessary 
in  order  to  avoid  the  difficulties  which  the  presence  of  a  square  root  introduces 
into  the  solutions  of  the  various  differential  equations. 

The  following  briefly  outlines  the  course  of  the  derivations  in  the  individual 
sections.  While  the  derivations  are  specifically  for  configuration  no.  1  the  re¬ 
sults  are  always  applicable  to  both  configurations  as  long  as  the  appropriate  fuze 

body  angles  are  used.  ’* 

i 

•( 

i 

.i 

GEOMETRY  OF  FUZE  BODY  CONFIGURATIONS  i 


The  following  section  contains  derivations  of  expressions  for  the  various  fuze  > 

body  angles  associated  with  the  pivot  holes  of  both  configurations  of  the  three  pass 
step-up  gear  mechanisms.  They  represent  extensions  of  previous  work  pertaining  to 
configuration  no.  1  given  in  reference  A-2,  appendix  A.  In  addition,  relationships 
between  the  unit  vectors  of  the  body-fixed  x-y  and  x'-y'  systems  (figures  A-l  and  A- 
2)  are  given  in  order  to  make  it  possible  to  derive  the.  escapement  relationships  in 
the  primed  coordinate  system,  whose  x'-axis  is  oriented  along  the  escapement  center- 
line  0  -0„  (ref  A-l).  Finally,  a  signum  function  is  introduced  so  that  common 

P  5 

programming  expressions  can  be  devised  for  both  configurations. 

Fuze  Body  Configuration  Ho.  I 

Figure  A-3  indicates  the  following  relationships  for  the  angles  ,  5^,  and 
for  fuze  body  configuration  no.  1. 

Angles 


From 


2  2  2 

(Rpj  +  r^)  ■  *  26?i(R2COSY2 


the  following  expression  is  obtained  for  ^ 


Y2  ”  cos 


-1 


(il  -  Up,  *  y) 

1  61 2 


Similarly,  from 


Y  3 


cos 


-1 


2 

<w2  + 


61, 


-  (R 


2  612  6*3 


El 


(A-l) 


(A-2) 
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Figure  A-3.  Fuze  body  configuration  no.  1  (three  pass') 


ie  following  is  obtained 


^3  =  ^2  +  y'3 


(A- 3  ) 


Lso  from 


Y4  =  cos 


«3  +  -  (pp3  +  rp4)2 


2  <?3  <U 


ie  following  is  obtained 

Y4  =  y3  +  Y4 


(A-4) 


(A- 5  ) 


ith  a,  the  distance  between  the  pallet  and  the  escape  wheel  pivots,  the  following 
s  obtained 


2  2 

4  +  <*'  5  “  2  4  6\  5cos  Y  5 


(A-6) 


hen. 


nd 


Y  5  =  cos 


6i  u  +  (Si , 


2  A  4  <R  5 


(A-7  ) 


Y5  =  Y4  +  Y5 


CA-8) 


Angles  6^ 


ince 


(R  ,  +  r  .) 

pi  p2 


+  -  2(Rpl  +  rp2)  6it  cos  6^ 


&2  =  COS 


-1 


(Rpl  +  rP2)  +  «!  -  «2 

2«1(Rpl  +  rp2) 


(A-9) 


imilarly , 


6a  =  cos 


-1 


(Rp2  *  cp3>2  *  ^ 

2  «2«p2  +  V 


(A-JO) 


36 


and 


the  following  is  obtained 


Angles 


With  equation  A-9 


( A- 1 1 ) 


(.A-12 ) 


(  A- 1 3 ) 


*  TT  -  d2  (A-14) 

Further,  with  equations  A-l  and  A-10 

e2  =  ^2  +  *  “  ( A- 1 5 ) 

With  equations  A-3  and  A-ll 

63  =  Y3  +  11  “  $4  (A-16) 

Finally,  this  furnishes  with  equations  A-5  and  A-13: 

84  =  Y4  +  «  -  65  ( A- 1 7 ) 

Figure  A-3  also  shows  that  the  angle  y'  between  the  x’-axis  and  the  unit  vector 
nj  is  given  by:  ** 

y"  m  x  —  e  ( A- 1 8 ) 
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where 


c  =  n  -  65  -  Y5 


(A-19) 


and  therefore 


y'  =  Sc  +  Ys  (A— 20) 

P  3 

The  unit  vector  n5  may  now  be  expressed  In  terms  of  the  primed  coordinate  system  as 
follows : 

n  c  =*  cos  y'i'  +  sin  f'j  '  1A-21) 

3  P  P 

Further,  the  unit  vectors  I'  and  j'  are  given  in  terms  of  the  x-y  system  by: 
i "  =*  005(84  -  it)i  +  sin( 84  ”  x)j 

i  "  35  -cos  84!  -  sinBiJ  (A-22) 

j'  =  k'  x  i'  3  -sin( 84  -  ir)i  +  003(84  ”  1T)  j 

j  '  =  sinB4i  -  cos  B4J  (A-23) 


Fuze  Body  Configuration  No.  2 

Figure  A-4  indicates  the  following  relationships  between  the  angles  yi , 
and  8j  for  fuze  body  configuration  no.  2. 


Angles  Y 

Since  these  angles  are  given  in  the  clockwise  direction  with  respect  to  the 
body-fixed  x  -  axis,  their  values  are  negative.  Thus,  with  the  help  of  equation  A- 
1: 


Y  2 


-1 

-cos 


+  «pl 

2 


+ 


(A-2A) 
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Y3  is  given  by  equation  A-2,  and  with  that  y  becomes 

3 

Y3  -  Y2  -  Y3  (A-25) 

The  negative  sign  for  Y3  is  necessary  in  order  to  make  y3  negative. 

Yi|  is  given  by  equation  A-4.  This  leads  to  the  following  expression  for  Y4  • 

Yi*  -  Y3  “  y£  (A-26) 

The  angle  Y5  is  again  taken  directly  from  equation  A-7 ,  and  with  that 

Y  becomes 
5 

Y5  =  Y4  -  Y5  (A-27) 

Angles  6^ 

The  angles  62,  63,  64,  and  65  may  be  taken  directly  from  equations  A-9,  A-10, 
A-ll,  and  A-13,  respectively. 

Angles  @1 

With  equation  A-9 

61  *  x  +  «2  (A-28) 

With  equations  A-10  and  A-24 

®2  =  Y2  +  n  +  53  (A-29) 

where  the  negative  sign  of  equation  A-24  for  Y2  is  accounted  for. 

Similarly,  with  a  negative  sign  of  y3»  one  obtains  with  the  help  of  equations 
A-ll  and  A-25: 

83  =  Y3  +  x  +  64  (A-30) 
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Finally,  with  equations  A- 13  and  A-24 


84  =  Y4  +  n  +  Ss 


(A-31 ) 


According  to  figure_A-4,  the  angle  y'  between  the  positive  x'-axis  and  the 
direction  of  unit  vector  ns  is  now  given  by  p 


where 


it  +  e 


e  =  it  -  6S  -  y's 


(A-32) 


(A-33) 


so  that 


Yp  =  h  ~  <$5  "  Y5  (A-34 ) 

The  unit  vector  n5,  when  expressed  in  the  primed  coordinate  system,  is  again 
given  by  equation  A-21,  i.e., 


ns  =  cosy^i'  +  siny'j  "  (A-35) 

Further,  the  unit  vectors  i“  and  j'  are  given  in  terras  of  the  x-y  system  by 
1'  =  cos(64  +  lr)i  +  sin(  84  +  ir)j 

i"  ■  -cos  64!  ~  sin  84J  ( A-36 ) 

j  '  *  k'  x  r  =  -sin(  84  +  tt)I  +  008(84  +  tt)J 

j"  =  81084!  ~  cos  84I  (A-37) 
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Common  Computational  Expressions  for  Both  Configurations 


In  order  to  find  common  programming  expressions  for  the  angles  of  both  config 
urations  (with  the  exception  of  the  angles  y'  of  equations  A-20  and  A-34),  th 
following  signum  function  s^  is  introduced,  wheife 

Sg  =  +1  for  configuration  no.  1  (A-38) 

=  -1  for  configuration  no.  2  (A-39) 

This  leads  to  the  following  expressions: 
according  to  equations  A-l  and  A-24 


-1 

Y2  =  s  6  COS 


«1  +  «2  -  <Rpl  + 
2  <R  j  d?  2 


(A-40) 


according  to  equation  A-2 

2  2 


-1 


«2  +  *I3  -  (a„2  +  rp3) 
2  2 


Y3  *  cos 

according  to  equations  A-3  and  A-25 


(A-41) 


T3  *  T2  +  s  6^3 


(A-42) 


according  to  equation  A-4 


according  to  equations  A-8  and  A-27 
Y5  =  Y4  +  s  6Y  5 

according  to  equations  A-9,  A-10,  A-ll,  and  A-13 


6?  =  cos 


-1 


(Rpl  +  CP2)2+  ~  A  2 

2  +  rp2^ 


5i  a  cos 


-1 


(R,2  +  ^  2  + 

2  «2«P2  +  rp3> 


6U  =  cos 


-1 


2  r,1  r>2 

+  fl3~ 

2  d?3<Rp3  +  rp4) 


fit  *  cos 


-1 


2  2  ^2 
L  +  <R  4  “  5 

2a  <(U 


(A-46) 


(A-47) 


(A-48) 


(A-49) 


(A-50) 


according  to  equations  A- 14  and  A-28 
-  it  -  s662 

according  to  equations  A- 15  and  A-29 
B2  -  Y2  +  *  ~  8653 

according  to  equations  A-16  and  A-30 
83  ■  Y3  +  *  -  a  664 

according  to  equations  A- 17  and  A-31 
B4  -  Y4  +  H  -  s  665 


(A-5 1 ) 


CA-52) 


(A-53) 


(A-54 ) 
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The  signura  function  will  also  be  useful  in  conditional  statements  which 
distinguish  between  the  of  configuration  no.  1  and  2  as  given  by  equations  A-2U 
and  A-34.  ^ 

Finally,  according  to  equations  A-22,  A-23,  A-36,  and  A-37,  identical  expres¬ 
sions  may  be  used  for  both  configurations  when  expressing  the  unit  vectors  of  the 
primed  coordinate  systems  in  terras  of  the  unprimed  ones,  i.e., 


1'  =  -cos  64!  “  sinfSiJ 


(A-55) 


j  "  =  sing4i  -  cos  64] 


(A-5b) 


DYNAMICS  OF  THE  PALLET  IN  COUPLED  MOTION 


While  the  following  derivation  for  the  dynamics  of  the  pallet  in  coupled  motion 
is  illustrated  with  figures  relating  to  configuration  no.  1,  all  expressions  are 
also  applicable  to  configuration  no.  2  as  long  as  the  appropriate  sign  for  the 
signum  function  s&  (eq  A-38  and  A-39),  as  well  as  the  correct  form  of  the  angle  y' 
(eq  A-20  and  A-34)  are  used.  ** 

As  stated  earlier,  the  dynamic  analysis  of  the  pallet  is  most  conveniently 
performed  in  the  primed  coordinate  system.  The  coefficient  of  friction  at  the 
pallet-escape  wheel  interface  and  at  the  pallet  pivot  has  the  designation  pj  . 

With  the  spin  u>  of  the  fuze  body  constant,  the  acceleration  of  the  center  of 
mass  of  the  pallet  is  given  by  figure  A-5. 

2  _  ,2  _  _ 

Acp  -  -id  A{5n5  -  (u>  +  *p)  rcp  [cos(ip+  \pc>i'  +  sin(ip+  ipc>j'] 

+  iprcp  [-sin(  ip  +  <p£  )i  "*  +  cos(ip+  ^c)j']  (A-57) 


This  simple  change  in  angles  is  all  that  is  required  since  the  relative 
rotation  directions  of  the  components  are  identical  for  both  configurations, 
and  therefore,  force  and  moment  expressions  retain  the  same  general  form. 
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With  the  above,  and  according  to  the  free  body  diagram  of  figure  A-5,  substitution 
into  Newton's  law  gives 


n  -pisJ’n^  +  F'i'-UiScF'  i'  +  F '  j  +  w is  '  j  ' 
nn  1  ^n  t  xp  l  yp  yp  y  xpJ 

“  m  [-u)  <^>505  -  ( oj  +  \p)  r^pCcosCtJi  +  +  sin(  iji  +  ) j  ') 


+  iJrCp(-sin(  ij)  +  ipc )  i  '  +  cos(  4»  +  4»c>j  ')  ] 


(A-58) 


The  signum  functions  s^  and  assure  the  proper  directions  of  the  friction  forces 
on  the  escape  wheel  -  pallet  interface  as  well  as  on  the  pallet  shaft,  regardless  of 
the  direction  of  pallet  rotation.  Thus, 


(ref  A-l,  eq  B-l)  and 


(A-59) 


Sc  «■  — (A— 60) 

Ul 

The  unit  vectors  n^  and  n^  are  adapted  from  equations  A-l  and  A-2  (ref  1). 

The  moment  equation  of  the  pallet  must  be  written  with  respect  to  the  acceler¬ 
ated  pivot  0^  ,  i.e., 

MoP  “  -r0p  x  mpr CP  +  HoP  (A-61) 

where 

Mqp  *  sum  of  all  external  moments  with  respect  to  pivot  Op 
rop  “  absolute  acceleration  of  point  Op 

Hop  »  time  rate  of  change  of  angular  momentum  of  the  pallet  with  respect  to 
point  Op. 


r 


Since  the  angular  velocity  w  of  the  fuze  body  is  constant  equation  A-61  takes 
the  following  form: 

2  _  _  _ 

Mop  =  -(-w  <>'5nb)  x  mprCP(cos(  ip  +  t|c)i'  +  sin(  tpf  <Jc )  j  ')  +  ipip  (A-62) 

Appropriate  computations  and  substitution  of  all  moments,  according  to  figure 
A-S,  furnish  the  final  moment  equation  (the  moment  arms  A^ ,  Bj,  Cj,  and  Dj  of  refer¬ 
ence  A- l  are  now  primed) 

~  2 

pn(Di  +  Cjy^)  -  +  Fyp)  =  Ipip  -  mprCp  d?  sin(  Yp-ip-.pc)  IA-63) 


where 


l)  I  —  c  cos(<>  -  a  -  ip)  (A-64) 

C  i  =  -  [r  +  c  sin(  <p  -  a  -  \p)  ]  (A-63) 


As  originally  discussed  in  reference  A-2,  Fxp  and  Fyp  represent  conservatively 
evaluated  pivot  force  components  which  assure  that  the  pivot  friction  moments  are 
opposed  to  the  rotation  at  all  times.  The  following  illustrates  how  this  goal  may 
be  accomplished. 

The  pivot  friction  components  Fxp  and  Fyp  are  first  obtained  from  the  following 
component  expressions  of  equation  A-58.  Subsequently,  they  are  transformed  to 
become  F  and  F  ,  respectively. 

*p  yp 

-P  sin(ip-a)  -  UiS4P  cos(ip-a)  +  F"  -  viiSeF' 
n  1  ^  n  xp  1  J  yp 

2  .2 

=  nip  [-ai  ^scos y'  -  (u+  ip)  rCpcos(ip+  ipc)  -  iprcpsin(ip+  ipc)  ]  (A-66) 

and 


P  cos(4>  -  a)  -  u,s UP  sin(<p  -  a)  +  F'  +  UiScF' 
n  1  4  n  yp  1  3  xp 

2  .2 

-  nip  [-ai  (WgsinYp  -  ( u>  +  ip)  rcpsin(ip  +  >pc)  +  iprcpcos(ip+  ipc)  ]  (A-b7) 
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Simultaneous  solution  of  equations  A-66  and  A-67  furnishes 


~  2  .  •  2  •• 

Fyp  =  PnAl  +  “  A2  *  2"|^J  *  a3  ±  ^4 

~  2  ^  .  2  •• 
Pxp  =  PnA&  +  “  A6  *  2  7^7  ^7  *  *  A7  ±  ^8 


where 

A  1  - 


A  2  = 


A3  = 


A  4  = 


As  “ 


A  r  = 


2 

Ul  (s  4-s  c,)sin(  (jr-a)  -  (I+S4S5U1  )cos(ij!-a) 


2 

1  +  pj 


mp[  ^s(  sin  y'-  H  is  5COS  y')  +  rcp(sin(  ij*^)  -  p  jS  geos  (  )  ] 

2 

1  +  Mi 


mprcp[sin(^«Pc)  ~  PiS  gcosC  i|^i|)c)  ] 
2 

1  +  Ml 


mprcp  (cos^)  +  MiS*in<  *►*<,)] 

2 

1  +  Pl 

2 

Ui(s  4-8  5)cos(  i)r-a)  +  (1+s  4s  5pi)sin(  ijr-a) 


2 

1  +  Mi 


Up  [-  <»s(cosYp+pis  gsinYp)  -  rCp(cos(  Ijrf-I()c)+Pis  5Sin( 


The  factor  w  / 1 co |  is  introduced  in  place  of  w  in  order  to  make  sure  that  the 
quantity  is  positive  regardless  of  direction  of  spin.  (This  too,  assures  th^t  the 
friction  moments  oppose  rotation.  The  driving  moment  is  proportional  to  w  ,  and 
therefore  independent  of  spin  direction.)  To  make  the  final  decision  concerning  the 
signs  of  equations  A-68  and  A-69 ,  these  forces  are  substituted  into  moment  equation 
a-63,  and  the  influence  of  the  direction  of  rotation  on  each  of  the  resulting  moment 
computations  is  explored: 


PJD1  +  CiU!S4  -  PpU^CA!  +  A5)]  ±  ppPis5u)  (A2  +  A6) 

2  2 

*  2P ph l  S5  j-~j  KA3+  A7)  ±  P pU i s 54*  (A3  +  A7)  ±  ppulS5iKA4  +  A8) 

2 

=  Ip'l'  "  mprcp  sin(Yp-»p-i(<c)  (A- 78) 


With  s^  positive  for  positive  rotation,  and  vice  versa,  while  all  other  para¬ 
meters  are  positive  at  all  times,  the  following  moment  components  of  equation  A-7R 
must  have  negative  signs  during  positive  rotation: 


-pnPpMlSs(Al  +  A5) 

2 

~P  pl*l  s5a*  <A2  +  Ae) 

.2 

"PpPl  s 5 4*  ( a3  +  a7) 


(A-79) 

(A-80  ) 

( A— 8 1 ) 


The  sign  of  the  term  containing  is  negative  and  is  controlled  by  the  sign 
of  iji  .  Therefore,  the  signura  operator  s^  is  omitted  and  the  term  becomes 


2 

"2PPU1  jff  *  (A3  +  A?) 


(A-82 ) 


The  choice  of  sign  for  the  friction  moment  term  in  equation  A-78,  which  is 
proportional  to  the  angular  acceleration  ifi  ,  is  discussed  in  detail  in  appendix  F. 
This  work  results  in  the  computational  rules  of  equations  A-88  and  A-89,  which  deal 
with  the  sign  in  the  effective  moment  of  inertia  lpR.  (Note  that  the  signum  func¬ 
tion  S5  function  has  now  been  omitted.) 
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With  these  considerations  equation  A-78  then  becomes 


2  a)  •  *  ^ 

PnA18  -  “  A19  -  2  T^y  ♦  A20  -  ♦  A2 1 

=  XPR*  '  mprCP  ****  isin(Yp  "  *  "  V  1 


(A-83) 


A  la  =  D  1  +  C  lV*  Is  4  -  PpPls  5<A  1  +  A  5)  (A'  84) 

A19  =  Pp  P  Is  5^ A  2  +  At>)  (A-85) 

A  2  0  =  PpPl^Ai  +  A  7)  (A-8b  ) 

A  2  1  =  Pp  U  [S  8(  A  -j  +  A  7)  (  A~8  7  ) 

I  =  I  +  A 22.  when  ip  and  1})  have  identical  signs  (A-88) 

PR  P 

IpR  =  Ip  -  A22’A_Z>  whea  ♦  and  'I'  have  opposite  signs  (A-89) 

A22  =  PpRiUi,  +  A8)  (A-90) 


Equation  A-83  is  now  used  to  find  an  expression  for  the  contact  force  Pn.  This 
jression  will  later  make  it  possible  to  establish  a  single  differential  equation 
r  the  escapement  in  coupled  motion.  Thus, 

.2  J  .  2  2 

IpR4H-A2  1  ^  +2-J^|-A  20^(0  Aig-mprCp  <R5u>  sin(  Yp-i|r-i|»c) 

=  *  ( A- 9 1 ) 

A  la 


*  The  parameters  An  are  not  sequential  at  this  point.  Ag  to  A^y  may  be  found 
in  equations  A— 111  through  A-116  and  equations  A-118  through  A-120. 


*  Care  must  be  taken  that  Ip  -  A22  does  not  become  negative.  If  this  occurs, 
Ipp  must  be  set  equal  to  zero.  See  footnote  A-8  where  the  possible  value  of 
lpR  in  the  free  motion  equation  of  the  pallet  is  discussed. 
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As  in  reference  A-l  the  variables  <j>  and  are  now  expressed  in  terms 
of  <*>  and  $  the  escape  wheel  variables  (ref  A-l,  eq  A- 15  through  A- 18  and  B-6) 


tp  =  U  4>  +  V 


(A-92) 


iji  =  Z  $  (Z  is  not  defined  in  ref  A-l) 


(A-93 ) 


where 


Z  =  — - - -  [Q  +  sin(  <J>  -  at)  } 

c  cos  ip  L  S  1 


(A-94) 


Equation  A-91  then  becomes 


p„  ‘  TTT,  I1™"*  +  ^  +  ipkv)‘2  +  *lSi  A^» 


2  2 

+  u)  A19  -  m(>rcp  <W5oj  sin(  Yp-ilr-^c)  ] 


(A-95) 


DYNAMICS  OF  ESCAPE  WHEEL  IN  COUPLED  MOTION 
(ESCAPE  WHEEL  INCORPORATES  PINION  NO.  4) 


While  figure  A-6  shows  a  free  body  diagram  of  the  escape  wheel  and  pinion  no.  4 
of  configuration  no.  I,  the  following  derivations  hold  for  both  configurations  as 
long  as  the  appropriate  signura  function  s^  is  chosen.  The  pivot  forces  Fx^  and  Fy^ 
as  well  as  forces  F-^  and  are  now  defined  in  the  general  x-y  system.  This  makes 
it  necessary  to  give  the  unit  vectors  nt  and  nn  in  terras  of  the  x-y  system. 

From  equations  A-l  and  A-2,  reference  A-l  (and  using  x'-y'  coordinates) 


nt  “  cos($  -  a)  1'  +  8ln(4>  -  a)  j 


(A-9b) 


n  "  -sin($-a)i'+  cos(4>  -  a)  j  ' 


(A-97 ) 


~ee 


Equations  A-55  and  A-56,  which  hold  regardless  of  configurations,  are  used  for 
the  following  transformed  expressions: 


nt  =  -cos(  $  -  a  +  84)  i  -  sin(  <ji  -  a  +  84)  j 


(A-98) 


and 


n 

n 


sin(  <j>  -  a  +  84)  i  -  cos(<j>  -  a  +  84)  j 


(A-99) 


The  expressions  for  the  unit  vectors  034  and  ^34  ,  as  used  in  the  analysis  of 
pinion  no.  4  in  reference  A-?,  section  A-la,  are  of  further  interest 


034  =  sin(03  +  83)  i  -  cos(83  +  63)  j 


(A-100) 


and 


°N34  =  COS  +  ® 3^  *  +  sin(B3  +  83)  j  (A-101) 

Using  the  concept  of  the  D'Alembert  force 

2 


T4  »  m4  <R4  id  (cos y 4  i  +  siny4  j),  (A-102) 

the  force  equation  of  the  escape  wheel  assembly  becomes  for  counterclockwise  rota¬ 
tion^  ^  (ref  A-2,  eq  A-35) 


_Pn  nn  +  MlS  ^n11!  +  F  34*  34  +  us  3F  34nN34  +  T4  +  F^i  +  wFy4i 

+  **x4I  '  FyJ  =  °  (A~103) 

Note  that  the  coefficient  of  friction  m  is  now  used  for  all  pivots  and  gear  tooth 
contacts  of  the  remainder  of  the  mechanism  train. 


A  3  See  appendix  F  for  description  of  motion  reversal,  i.e.,  clockwise  escape 
wheel  rotation.  This  may  occur  after  severe  impacts. 


A-6 


The  signum  functions  Sp  82,  and  83  are  defined  in  reference  A-2  in 
connection  with  the  tooth  contact  friction  of  the  various  meshes. 
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The  moment  equation  of  the  escape  wheel  for  counterclockwise  rotation  is  writ 
ten  with  the  help  of  figure  A-6 


"Pn^Al  +  +  Py4>  +  rb4F  34  "  ^3^d3  “  a  3^  F  34  “  14<t»  (A-104) 

where 

At  “  b  cosa  +  g  (A-105) 

> 

B]  *  b  sina  (A-106) 

The  pivot  forces  F  and  F  are  treated  again  in  the  manner  discussed  ear 
lier.  They  are  obtained  rrom  th?  component  expressions  of  equation  A-103,  i.e., 

-P^sinC  <f>—  art- Bi*)  -  s  4y  jP^cosC  <^ert-84)  +  F  34sin(  8 3+63) 

+  ns  3F  34cos(  8 3+0  j)  +  T4cosy4  +  F^  +  gFy4  =  0  (A-107) 


and 


P^cosC  <t>-ori-84)  -  s  4pjPnsin(  <}r—  ort- 34)  -  F  34cos(  83+93) 


+  os  3F  34sin(  83+63)  +  T^lny^-  Fy4  +  uF^  *  0  (A-108) 

Simultaneous  solution  of  equations  A-107  and  A-108  gives 

Fy4  -  PnA*+F3,A10+T4Au  (A-109) 

fx4  *  pn  A  1 Z*"  F  3  4  A  1  3  +  T4  Al4  (A-110) 
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Aq  = 


-(  M]S  4~p)sin(  4>—  or*- p 4)  +  (  1+UU  18  q)C08(  I^ori-Bq) 


1  +  U 


A10 


-u(l-s  3)sin(  $31-03)  -  (1+u  s  3)cos(  $3f  03) 


1  +  n 


An 


sitiY4  "  ijcos y 4 


1  +  u 


12 


(  1+UlUS  4)sin(  <(rcd-$4)  +  (s  4lJ  pli)cOs(  (|r-oH-$4) 


1  +  u 


A 1  3 


-(l+u  s  3)sin(  $3+-03)  +  u(l"S  3)cos(  $3+03) 


1  +  U 


A  1 4  “ 


-C0SY4  **  usinY4 


1  +  M 


(A-lll) 


( A- 1 1 2 ) 


( A- 1 1 3  ) 


(A-114) 


( A- 1 1 5 ) 


C  A— 116) 


Substitution  of  the  pivot  force  component  expressions,  with  the  indicated 
signs,  into  the  moment  equation  A-104  will  furnish  the  correct  pivot  friction  mo¬ 
ments.  The  resulting  expression  is  then  solved  for  the  contact  force  P 


-1 4  (j>  +  F34  A15  -  T4  Ai6 


p  - - - - 

n  A17 

( A— 1 1 7  > 

where 

A 15  *  rb4  "  y  [83<d3  “  *3)  +  P**(A10  +  A13)] 

(A-118) 

A  16  ”  UP  4(  A  1 1  +  A  14) 

( A- 1 1 V  ) 

A17  -  A  j  +  U184BJ+  UP  4(  A  g  +  A12) 

(A-12U) 

COMBINED  COUPLED  MOTION  DIFFERENTIAL  EQUATION 
FOR  ESCAPE  WHEEL  AND  PALLET 

(APPLICABLE  TO  BOTH  CONFIGURATIONS  WITH  APPROPRIATE  s6) 


To  obtain  a  combined  coupled  motion  differential  equation  for  the  pallet  and 
escape  wheel  In  terms  of  the  escape  wheel  angle  (J>,  equations  A-95  and  A-117,  both 
for  the  contact  force  Pn  are  set  equal  to  each  other.  This  furnishes 

2  2  2 

|  f  i+A  i  8  +  IpRA  i  7U  J  <j>  +  j^IpRA17V  +  A17A21Z  j  +  2-j^y  A17A20Z(£ 

2 

■  A  1  &A  18*"  34  “  A16A18T4~  A17o)  j^A  ^9  -  roprCp  5sin(  J  (A- 12 1 ) 

To  complete  the  system  differential  equation,  it  Is  necessary  to  find  an  ex¬ 
pression  for  the  contact  force  F^  in  the  above.  This  will  be  accomplished  by 
combining  the  appropriate  differential  equations  for  gear  and  pinion  no.  2  and  3  and 
for  the  rotor. 


DYNAMICS  OF  ROTOR  (GEAR  NO.  1) 
(APPLICABLE  TO  BOTH  CONFIGURATIONS  WITH  PROPER  s6> 


A  free  body  diagram  of  the  rotor  for  configuration  no.  I  Is  shown  in  figure  A- 
7.  The  acceleration  A  of  its  center  of  mass,  with  constant  spin  velocity  u>  ,  Is 
given  by 


2  ,  2  [-  _  _-i 

A„„  •  “to  ^1  -  (urt-^j)  rcl  |cos(  tlc+<h)i  +  sln( <t>lc+* i> j  J 


CR 


+  <{.1rcl  -sin(  <t>lc+4>1)l  +  cos(-}>lr+ 


'1C 


>1>5  ] 


(  A- 1 2  2  ) 


Since  the  motion  of  the  rotor  must  be  expressed  in  terms  of  the  escape  wheel 
variables,  the  following  transformations  are  now  introduced: 


where 


■  N(,  j  l 

(A- 123) 

4> i  ”  4> 

(A-124) 

NP4  NP3 

Np2 

N,tl  N  N 

N„. 

(A- 1 23 ) 
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The  rotor  angle  4>.  +  44  ,  Is  expressed  as  follows: 


♦lC  +  *1  “  *1C  +  N41  *T 


(A- 12b) 


where  4>T  represents  the  total  rotation  of  the  escape  wheel  from  the  inception  of 
motion.  (The  section  on  Geometry  of  Fuze  Body  Configuration  describes  the  manner  in 
which  is  obtained  as  a  function  of  the  instantaneous  angle  $.) 

Equation  A-122  then  becomes 


2  .2 

krf)  *  -w  (S^i  -  (wbNm4>)  rr)  [cos(4>lc+N41(j)T)i  +  sin(  <fr)r+N4  i<|>T)j 


1 C  4 


+  Nm*  rcl  [-sin(*1(,bN41<Jir)I  +  cos(<Jilc+N414)T)j  ] 


(  A- 1 2  7  ) 


Newton's  force  equation  is  now  written  for  clockwise  rotation  of  the  rotor  with 
the  help  of  equation  A-127  and  the  free  body  diagram  of  figure  A-7:^-^ 


-Fi2nI2  -  u  slFl2%i2  ~  Fxli  +  bFylI  +  FylJ  +  \*xli  -  ®  iAc 


(A-128) 


where 


ni2  =  sin(Bi  +  0 i)i  -  cos(8!  +  G^j 


(A-129) 


n  .  =  cos (  B 1  +  0 i>i  +  sinCBi  +  G^j 


( A- 1 30  ) 


(ref  A-2,  eq  A-78  and  A-79) 
Further,  in  figure  A- 7 


F21  *  _F  12°  12 


(A-131) 


Ff9i  “  -  M  Sl  F12  n 


( A- 1 3  2 ) 


(ref  A-2,  eq  A-103  and  A- 104) 


Description  of  motion  reversal,  appendix  F. 
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The  moment  equation  must  be  written  in  the  manner  of  equation  A-61,  with  re¬ 
spect  to  the  accelerated  point  0^.  The  pivot  reactions  and  F  ^  are  treated  so 
that  the  associated  friction  moments  retard  the  clockwise  rotation  of  the  rotor. 
This  leads  to 


RblF12R  “  US  xa  XF  12k  +  up  i(Fxl+Fyl)k 

2 

=  -(-u)  (SJji)  x  ra  jrcl  [cos(^1G+N41<j>T)i  +  sin(  j  ] 

+  I  iNm.jk  ( A- 1 3 3 ) 

This  becomes 


RblF  12  "  ^  la  1F  12  +  bPi<Fxl  +  Fyl> 

2  , 

»  m  jo)  !rclsin(  i()lc  +  N,,^)  +  I XN 41  <p  (A-134) 

The  forces  F  ^  and  F  are  obtained  from  the  simultaneous  solution  of  the 
following  component  expressions  of  equation  A-128: 


-F  12sin(  -  u  s  jF  12cos(  B  i+0  i)  -  F^  +  pFyl 

2.2 

=  mx  [-u>  -  (orFN41^)  rclcos(  <(>lc+N41ij)T)  -  ]  (A-135) 

F  12cos(  8  j+6 j)  -  us  iF12sin(  Bi+Qj)  +  F  +  uF^ 

,  .  2 

3  mi  [-(urt-N41<|>)  rclsin(  ^C+Nm^)  +  N4i(Jirclcos(  4>1C+N4l  ]  (A-13b) 


Simultaneous  solution  of  the  above  gives 

2  2  .2 

Fyl  “  F  12^23  ±  u  A24  ±  -1^4 1^25  ±  <N41$)  A25  *  N  4 1^20 


and 

2  2  2 

Fxj  “  F12A27*  '0A28±  2-j~N  m  2  9  *  <N41$>  A29±  N41<a3U 


(A-137) 


(A- 138) 
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where 


A 23  * 

A 24  * 

A  2  5  = 

A  26  = 

and 

A  2  /  = 

a28  = 

A28  = 

A  in  ■ 


2 

p(l+s  i)sln(  B  i)  -  (1-p  s  i)cos(  6 1"4-  0 1 ) 

2 

1  +  P 

n»i[u  +  rcl(  ycosC^^+Nm^J+sinC^^+N^j^))  ] 

2 

1  +  P 

ra  1r(,1  [pcos(  <f>lc+N4l^r)+sin(  <i>lc+N4  lij.)  ] 

2 

l  +  p 

mird  [cos(^c+N4l<hr)-lBin(4.lc+N41<Hr)  ] 

2 

1  +  p 

2 

( 1-p  s  !)sin(  B  i+0  !>  +  p(  1+s  i)cos(  B  i+e  i) 

2 

1  +  p 


nll[<1{l+  L  C  cos  C  d>l  C+N  4  !  4^)  -  psin(  <tilc+N4l({»T))  ] 

2 

1  +  P 


mircJcosC  l>ic+N^V  "  lBln(  *ic+N,4l<V  I 


1  +  P 

mlrcl  [sln^'tllC+N4  1<hr)  +  PCOs(  lJ'lc+N4l\)  ] 


(A- 139) 


(A- 140) 


( A- 1 4  1 ) 


(A- 142) 


(A-143) 


(A-144) 


(.  A- 145 ) 


(.  A—  1 4(> ) 


1 


Equations  A-137  and  A-138  are  now  substituted  into  the  moment  equation  A-134  as 
follows : 


RbjF  12  “  P8  la  1F  12  +  UP  l[F  12(a23+A27)  *  U(A24+A28) 


.  2 


±  2-t~t  N  4  i|(A  25^A  29)  *  ( N  4 1  $)  (A25FA29)  ±  Nill<t>(A2e+A3o)  ] 


m iio  d^r  sin(  <t>lc+N4l<hr)  +  1  lN 4 1 4> 


(A-147) 


In  order  to  have  the  pivot  friction  moments  positive  for  clockwise  rotation, 
the  following  terms  must  be  positive: 


U  Pi  F  (A2j  +  A 2 7) 


2 

M  Pjiv  (A24  +  A2fl) 


U  Pi 


.  2 

(N4i<j>)  (A25  +  A2y) 


The  term 


2 

2  M  pl  N4l^A25fA29^ 


must  also  be  positive  for  clockwise  rotation  of  the  rotor.  This  can  only  be  accom¬ 
plished  by  making  the  term  negative  since  is  negative.  To  avoid  the  use  of  p  as 
a  signum  function,  as  required  for  motion  reversal  according  to  appendix  E,  its 
value  is  made  absolute.  Possible  motion  reversal  is  then  only  affected  by  the  sign 
of  $  ,  the  angular  velocity  of  the  escape  wheel,  which  is  normally  positive.  Then 
the  resulting  term  has  the  following  form: 


2 

"2  ^A25  +  a29) 

In  order  to  determine  the  sign  of  the  pivot  friction  moment  which  is  propor¬ 
tional  to  the  angular  acceleration  $  of  the  escape  wheel  in  equation  A-147,  the 
ideas  presented  in  appendix  F  are  used.  To  be  able  to  do  this  let  the  value  of  the 
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coefficient  of  friction  of  this  term  become  absolute,  so  that  It  ceases  to  serve  as 
a  directional  signum  function  In  the  sense  of  appendix  E.  Further  let  the  expres¬ 
sion  be  changed,  for  the  time  being,  so  that  it  becomes  a  function  of  the  rotor 
acceleration  ^  .  With  the  above  and  In  the  sense  of  equation  F-2,  appendix  F,  the 
absolute  value  of  this  friction  moment  may  be  expressed  as 

MAA  “  |u|Pi  (A26  +  A30)  4> l 

From  this  point  on  one  may  use  the  reasoning  of  appendix  F  directly,  keeping  in 
mind  that  |u|Pi(A26  +  A30)  and  <J> 3  are  now  used  instead  of  A 22  and  i|>,  respectively. 

As  for  the  four  cases  in  appendix  F,  the  effective  moment  of  inertia  IjK  takes 
two  forms: 

I  =  Ii  +  |  u  |  P 1  (A26  +  A30),  when  <(> !  and  41 1  have  the  same  sign, 

IK 

and 

I1R  =  1  1  ~  I  Pl  (A26  +  A  30),  when  <)>i  and  <j>j  have  opposite  signs 


Equation  A-147  gives  all  relevant  expressions  in  terms  of  the.  escape  wheel  var¬ 
iables  ifi  and  .  Since  they  both  are  proportional  to  <f)j  and  by  the  identical 
gear  ratio  (or  in  the  case  of  the  two  pass  train  by  way  of  N^j),  one  may  readily 
extend  the  above  computational  rule  to  the  escape  wheel  variables.  This  is  re¬ 
flected  by  equations  A- 155  and  A- 156. 

The  above  considerations  give  the  moment  equation  A-147  the  following  form: 


F12  [Rbl  "  US  ia  1  +  UP  i(A  23^A  27)  ]  +  w  P  lb(A  24"*"A  28> 

J-  .2 

“  2  |  p  |  P  1  -j~~  N 1*  i  A  2  5*"A  2  9)  +  UP  i(  N  m  <t>)  (A25+-A29) 

2 

=  mltcl  ^l«o  sin(  4>1C+Nm<fcr)  +  C I  !  +  UP  i(A  26+a3U))N  4  ^  (A-148) 
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Finally  the  above  expression  is  solved  for 


.  .  2 

-A3/FN4lA33|  -  A3tl(Nltl|)  +  A3^in(  $ .  +Ni)1i( l,)  +  ^ i  KN *+  i 4> 

K - - - - - 1—1* -  (  A- 1 4  9  ) 

Ail 

wliere 


A  3  i  =  R,.  -  VB  ja  t  +  UP  i(A  23  +  A27)  (A-150) 

2 

A 32  =  id  P1UIA24  +  A2a)  C A— 151) 


2 

a33  =  2  |P|bl  <A25  +  A2 3) 


(A-152) 


A  3U  "  PP  1  ^A  25  +  A  2y) 


( A- 1 5  3 ) 


2 

A  3  5  »  m  l  ^  irc  L  a) 


(A-154) 


l^R  =  lj  +  |  ul  P  1  (A  26  +  A  30) ,  when  |  and  $  have  the  same  signs  (A-155) 


I^R  =  Ij  -  iuIpi  (A26+  A  30)  >  when  |  and  4>  have  opposite  signs  (A-156) 


DYNAMICS  OF  GEAR  AND  PINION  SET  NO.  3 
(APPLICABLE  TO  BOTH  CONFIGURATIONS  WITH  PROPER  s&) 

The  moment  equation  of  gear  and  pinion  no.  3  may  be  taken  from  reference  A-2  if 
its  right  hand  side  is  set  equal  to  1 3413  rather  than  to  zero.  (The  dynamic  force 
equation  was  already  satisfied  by  the  inclusion  of  the  inertia  force  T3.)  Thus, 
with  the  help  of  equation  A-65,  reference  A-2,  and  letting 

<t> 3  -  N43  $  (A- 157) 
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( A— 1 58) 


i 


N43 


N 

N 


P4 

G3 


following  expression  is  obtained 


3F34  -  MS  3a  3F  -  rb3F23  +  MS2(d2  “  a2^F23 

+  bP  3  [F  2  3( A  3$  +  A 3 9^  +  T3<A37  +  A4<))  +  F34<A38  +  A4l>  ] 


~  I3N4341 


(1  +  M  s2)cos(02  -  02)  +  m(s2  -  l)sin(  02  -  02) 


1  +  u 


sin'y’3  +  ucos  Y3 


1  +  u 


( 1  —  u  s3)cos(03  +  03)  -  uC  1  +  S3)sin(  63  +  03) 


1  +  M 


(1  +  M  s2)sin(02  -  02)  +  u(l  -  s2)cos(02  -  02) 


1  +  M 


UsinY3  -  cos y 3 


1  +  M 


2 


(L  -  M  83)310(63  +  03)  +  m(1  +  s3)cos(03  +  O3) 

2 

1  +  U 


(A- 159) 


(A-160) 


(A-161) 


(A-lb2) 


(A-163) 


(A-lb4) 


(A-165) 
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[f  equation  A-159  is  solved  for  ^23 »  c*ie  following  is  obtained: 
FJ‘*A42+T3A43  -  I  3  N4  3 

F  2  3  = - - a -  CA-I66) 

A44 

where 

A42=Rj,3_lJ[s3a3-p3^A3B+A4l)]  (A- 167) 


A4  j  =  UP3(A3/  +  A  40)  (A-168) 

A  4  4  =  rb3  -  M[s  2(d  2  "  a  2)  +  P3^a36  +  A  iy)  ]  (A-169) 


DYNAMICS  OF  GEAR  AND  PINION  SET  NO.  2 
(APPLICATION  TO  BOTH  CONFIGURATIONS  WITH  PROPER  s&) 


The  moment  equation  of  gear  and  pinion  set  no.  2  may  also  be  taken  from  refer¬ 
ence  A-2  with  the  provision  that  its  right  hand  side  is  set  equal  to  l2t|>,  rather 
than  equal  to  zero.  (Again,  the  force  equation  was  already  satisfied  by  the  inertia 
force  T2.) 

With  the  help  of  equation  A-92,  reference  A-2,  and  letting 
$2  =  ^42  ( A— 170) 


where 


N42  = 


P4 


G3 


N 

N 


P3 

G2 


the  following  is  obtained: 


(A-171) 


"Rb2F23+  Ps2a2F23+  rb2F  12  _  1,8  l^d  l_a  1^F  12 


-  UP2  (F  12^a45FA4q)  +  T2(A4gfA49)  +  F  23(A4 yfA 5y)  ] 


1 2N  42$ 


(A-172) 
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where 


u(  1  -  s  i)sin(  8 1  +  9 1)  +  (1  +  M  i)cos(  8 1  +  0  i) 


1  +  u 


sin Y2  "  UCOSYz| 


l  +  u 


u(  1  +S2)sin(82  -  62)  +  (1  ~  Vs  2)cos*-  02  —  0  j) 


A47  = 


1  +  U 


u(  1  -  s  i)cos(  81+  @i)  -  (1  +  us  1)  sin  (81+  0 1) 


A48  = 


A49  = 


1  +  u 


UsinY2  +  COS  Y2  I 


1  +  u 


(1  -  U  s2)sin(82  “  ®2)  “  hd  +  s2)cos(82  ”  62) 


1  +  u 


The  foliowing  is  found  if  equation  A-172  is  solved  for  F^: 


F  23  A  5 1  +  T  2  A  52  +  I  2  N4  2 


where 


A51  =  H,  ,  “  U  [s  232  _  P 2^A 47  +  A5u) 


A  52  =  UP 2^A 45  +  A49) 


A53  =  ru,  “  ^  (  3  l<d  l  "  a  l)  +  P/A45  +  A4b>  ] 


(A-173) 


(A- 174) 


(A-173) 


(A-176) 


( A- 1 7  7 ) 


(A- 178) 


(A-179) 


(A-180) 


C  A— 181) 


(A- 182 ) 
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DYNAMICS  OF  COMBINED  SYSTEM  IN  COUPLED  MOTION 
(APPLICABLE  TO  BOTH  CONFIGURATIONS  WITH  PROPER  s&) 

In  order  to  obtain  a  single  differential  equation  for  the  total  system  in 
coupled  motion,  in  terms  of  the  escape  wheel  variable  <(>,  an  appropriate  expression 
for  the  contact  force  F^  must  be  substituted  into  equation  A— 121.  To  this  end, 
equations  A-149  and  A-179,  both  expressions  for  the  force  and  originating  from 

the  expressions  for  the  rotor  and  gear  and  pinion  set  no.  2,  are  first  set  equal  to 
each  other  and  then  are  solved  for  the  force  F23: 

^23  =  A)^AJ7  t  ~  A  3  ;A  5  2T  2  +  A  35A53sin(  4>lc+N4l4>T) 

.  2  .2 

+  A  3  jA  5 3N  i*  3 <p  -  A  34A  53N  43<j» 


+  (A53I1rN41  -  A  31I  2N  42>  4>]  (A-I83) 


The  just  determined  expression  for  F23  is  set  equal  to  equation  A-166  for  , 
which  originates  from  the  dynamics  of  gear  and  pinion  set  no.  3.  This  results  in: 


F  34  =  A  3iA42A5!  -A32A4<+A5 3"  A31A44A52T2-  A31A43A51T3 


+  A  3sA44A  53Sin(  i})1c+N4  jij^)  +  A  3  3A  44As3N4  t  <f> 


.  2 

-  A  34A44A  53(N4  jj)) 


+  (A44A  53IjrN4  1  -  A  3  iA44I  42  +  A  3  iA &11  3N  4 3)  ]  (A-184) 
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Finally,  the  above  is  substituted  Into  equation  A-121.  The 
differential  equation  for  coupled  motion  is  obtained: 


••  .Z 

A54<t>  +  A  S5<{> 


+  A56$  =  A57  +  A58sin(  yp-tjr-il>c)  +  A5ysin(  ^  C+N4 


where 


A  I5A  1^4^53 

a54  =  l*Aia+  1prai^  -  ■a31a42a51'  N4iI1R 

A15Ala 

2  A  15A  18A34A44A53  2 

Ass  =*  iprAi7v  +  a17a21z  +  aT^aH  Ni*i 

w2  A15A  ihAjjA^iAsj 

Asb  =  2— A17A2oZ  -  -T~A~7ir~  Nl41 

A15A18A43  A15A18A44A52 

As?  =  “A  18*18*4  T3  a^aJ7 


A  1  5A  18A  32A44A  53  2 

- - r — 7 - 7 -  "A  17A  13O) 

A  3  lA  4  2A  5 1 


A  5b 


A  17“  mprCF  <?5 


Asy 


A  isA  19A  35A  44A  53 


following  system 

(A-185) 

( A-186 ) 

CA-187) 

(A- 188) 

(A-189) 

(A-19U) 

(.  A— 191  ) 


A  3  1A  4  ^A  5  1 


CONTACT  FORCE  EXPRESSIONS  FOR  COUPLED  MOTION 
(APPLICABLE  TO  BOTH  CONFIGURATIONS  WITH  PROPER  s6) 


The  following  force  expressions  for  various  contact  forces  are  useful  for 
strength  calculations,  and  are  applicable  to  both  configurations. 

According  to  equation  A-184 


F  34  " 


1 


A  3  ^4^51 


[  “^32^44^53 


A31A44A52T  2 


A31A43A51T  3 


+  A35A44A  53Pln(  <t'ic+N41(kr)  +  A  3  3A  44A  S3N  4 1  + 


.  2 

“  A3itA44A53(N41<(>)  + 


(A44A53I1RN41 


-  A31Am(I  ;J<  42  +  3N43)*t> 


iA-192) 


According  to  equation  A-166 


F34A42+T3A43~  I3N4  3<i 


23 


L44 


According  to  equation  A-123 


12 


F2jA  51  +  T2a52+  1  2^*42  $ 


■53 


(A- 193) 


(A- 194) 


The  contact  force  P  ,  between  escape  wheel  and  pallet,  may  either  be  expressed 
in  terras  of  the  escape  wneel  variable  $  according  to  equation  A-117,  or  in  terras  of 
the  pallet  variable  ^according  to  equation  A-91.  Thus, 
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(A- 195) 


-14<f>  +  ~  T4a16 


‘17 


or 


Pn^  =  I^  [  *PR*+  a21^+  2-jfj  A20|+  <A19 


-  mprCp  ^5“  8in(Yp-<|r-*c) 


(A-196) 


DIFFERENTIAL  EQUATIONS  FOR  FREE  MOTION  REGIME 
(APPLICABLE  TO  BOTH  CONFIGURATIONS  WITH  PROPER  s6) 


Both  the  differential  equations  for  the  free  motion  of  the  pallet  and  the 
system  composed  of  the  escape  wheel,  gear  and  pinion  no.  3,  gear  and  pinion  no.  2, 
and  the  rotor  may  be  obtained  from  certain  coupled  motion  expressions. 


Free  Motion  of  Pallet 

Let  Pn  =  0  in  equation  A-91.  After  some  rearrangement,  the  following 

differential  equation  results: 


.2 

Aeo'l'  +  Aa +  Afei’J’  =  "  a62  +  Aea8111^  ~  t  -  <J>C) 


i A-197 ) 


where 


A60 


(A-19B) 


For  free  motion,  IpR  cannot  be  zero  since  it  would  make  the  value  of  in 
the  Runge-Kutta  solution  indefinite  (footnote  A-A ) . 
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Acu-  2T:7A20 


(A- l 99) 


A  62  "  w  A  i  vj 


(A-200) 


A  63  “  W  mpvcV  <l5 


Free  Motion  of  the  Escape  Wheel,  Gear  Train,  Rotor  System 


(A-201 ) 


For  the  present  case  start  by  letting  P  -  0  in  the  escape  wheel  expression  A- 
117,  i.e.,  . . . 


1 4 ^  “  Ai  5F  34  -  T4A16 


(A-202) 


Equation  A- 184 ,  which  furnishes  an  expression  for  force  F34  ,  is  then  substituted 
into  equation  A-202.  This  provides  the  desired  second  free  motion  differential 

equation 


.2 

A g, 4$  +  A g 5<j)  -  A 6 &4>  -  A67  +  AggsinC  i})^  +  N 4  1 4>rj, ) 


(A-203) 


whe  re 


A  64 


I  4 


AlsA4iAs3  A15A44. 

A31A42A51  N,)lIlR  +  a42a51  n^12  “ 


A  1 5 
A42 


N43I3 


(A-204 ) 


A  65 


A15A34A44A53  2 

A  3  1A  4  'A  5  1  N  4 


(A-205) 


A  6  6 


A15A33A41A53 

A  3  lA  4  ^A  5  1 


(A-20b) 
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(A-207  ) 


A15A43  A15A44A52  a  15^  31^44^53 

a67  =  -AI6T4  "  A42  T3  -  Al42A51  T2  "  Aa^As;  " 


AlsA  3^44^3 
168  a31a42a51 


(A-208) 


CONTACT  FORCE  EXPRESSIONS  FOR  FREE  MOTION 
(APPLICABLE  TO  BOTH  CONFIGURATIONS  WITH  PROPER  s&) 

The  contact  force  Fp^4  may  be  obtained  from  equation  A— 117  by  setting  Pn  equal 
to  zero,  so  that 


I  4<t>  +  T  t,A  16 
FF34  =  A ~S 


(A-209) 


The  expression  for  the  contact  force  Fp23  has  the  same  form  as  equation  A-193 
since  it  originates  from  the  dynamics  of  gear  and  pinion  set  no.  3 


FF34A42  +  T  ^**3  "  1  3N43<t> 


‘  F23 


l44 


(A-210) 


Similarly,  the  form  of  Fp^2  is  like  that  of  equation  A-194 
FF23A  5  l  +  T  2*  52  +  1  2N  424* 

FF12 


(  A-  2  1 1 ) 


CHANGES  IN  IMPACT  EXPRESSIONS 
(APPLICABLE  TO  BOTH  CONFIGURATIONS  WITH  PROPER  s6) 


While  generally  the  impact  description  of  reference  A- 1  remains  unchanged,  the 
total  moment  of  Inertia  Ig-jQT  of  the  escape  wheel  is  now  increased,  since  it  must 
also  account  for  the  presence  of  the  gear  train  and  the  rotor.  Therefore, 
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(A-212) 


2  2  2 

ISTOT  “  r*»  +  T3N43  +  l2N“2  +  IlN41 


where 


=  escape  wheel  -  pinion  no.  4  moment  of  inertia 
X^  =  Rear  and  pinion  set  no.  3  moment  of  inertia 

X^  =*  gear  and  pinion  set  no.  2  moment  of  inertia 

X^  =  rotor  moment  of  inertia 

See  equations  A-125,  A-158,  and  A-171  for  the  transmission  ratios. 
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APPENDIX  B 

DYNAMICS  OF  ROTOR  DRIVEN  S&A  DEVICE  WTTH  A  TWO-PASS 
INVOLUTE  GEAR  TRAIN  AND  A  PIN  PALLET  RUNAWAY  ESCAPEMENT 


i 


DESCRIPTION  OF  SYSTEMS  AND  OUTLINE  OF  DERIVATIONS 


Figures  B-i  and  B-2  show  the  two  configurations  of  the  S&A  mechanism  discussed 
in  this  appendix.  These  devices  differ  from  those  described  in  appendix  A  because 
they  contain  only  two  step-up  gear  meshes.  Since  the  escape  wheel  still  is  assumed 
to  have  counterclockwise  rotation,  the  motion  of  the  rotor  must  also  be  counter¬ 
clockwise,  and  the  gear  and  pinion  set  no.  2  must  turn  in  a  clockwise  direction. 
This  requires  certain  changes  in  the  rotor  geometry  as  well  as  in  the  direction  of 
the  lines  of  action  of  the  gear  contact  forces  F^  and  F2-J.  Naturally,  all  quanti¬ 
ties  previously  associated  with  gear  and  pinion  set  no.  3  are  now  absent. 

The  general  approach  to  the  modeling  of  these  devices  is  the  same  as  that  used 
in  appendix  A.  While  the  derivations  are  specifically  for  conf iguration  no.  1,  the 
results  are  applicable  to  both  configurations  as  long  as  the  appropriate  fuze  body 
angles  are  used. 


GEOMETRY  OF  FUZE  BODY  CONFIGURATIONS 


Fuze  Body  Configuration  No.  1 

Figure  B-3  indicates  the  following  relationships  for  the  angles  y  ,  <5^ 
and  for  configuration  no.  1: 

According  to  equation  A-l  (app  A) 


-1 

y  2  =  cos 


+ 


fl  , 


(R 


£i 


+  V 


2  (R  j  (f?2 


Further,  according  to  equation  A-2  (app  A) 


IB-1) 


Y3  =  cos 


-1 


«2+  <?3  -  (Rp2  +  rp3>2 


2  (R2 

Finally,  from  equation  A-3  (app  A) 


(B-2) 


Yi  =  Y2  +  Y3 


(B-3) 


FriEGSDLNG  PADS  BUNK -NOT  FILMED 


L 
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gure  B-3.  Fuze  body  configuration  no.  1  (two  pass) 


The  angle  y£  must  be  newly  derived.  Thus, 

2  2  2' 

_!  ^3  +  ^4  -  a 

yC  =  cos  - y-r— - 

L  A<  3  (1?  4 

Then,  as  In  equation  A-5  (app  A) 

1 4  =  Yj  +  YC 


(B-4) 


(B-5) 


Angles  6^ 

The  angles  6 2  and  63  can  be  taken  directly  from  equations  A-9  and  A-10 
(app  A),  respectively,  i.e., 


&)  =  cos 


-1 


(Rpl  +  cp2>  2  +  flI  "  *'2 


2  *4(Rpi  +  rp2^ 


and 


5 3  =  cos 


(Rp2  *  r,3>2+  *'2-  fl2 


2  «2  <Sp2  +  rp3> 


(8-6) 


(  B— 7 ) 


The  angle  64  must  take  the  pivot-to-plvot  distance  a  of  the  escapement  into 
account.  Therefore, 


64  =  cos 


2 

a  + 


2a  <R3 


(B-8) 


Angles 

The  angles  83  and  82  may  also  be  taken  directly  from  equations  A-L4  and  A- 
15  (app  A),  respectively,  i.e., 


B 1  =  it  -  (B-9) 

and 

02  =  Y2  +  "  <5 3  (B-10) 
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The  angle  83  is  found  with  the  help  of  angle  64  of  equation  B-8.  Then,  similar  to 
equation  A-16  (app  A) 


83  =  Y3  +  ^  -  64 


(B-ll) 


Figure  _B-3  shows  the  angle  y'  between  the  positive  X-  axis  and  the  unit 
vector  n^.  It  is  given  by  P 


y"  =  IT  -  £ 

p 


where  e  is  obtained  with  the  help  of  equations  B-4  and  B-8 


(B-12) 


£  =  ti  -  54  -  YC, 


(B-13) 


and  therefore 


y'  -  64  +  YC  (B-14) 

P 

The  unit  vector  n4  is  expressed  in  terms  of  the  primed  coordinate  system  as  follows: 

n4  =  cosY^i^  +  siny'j  ■*  (B-15) 

P  P 

Further,  the  unit  vectors  i "  and  j  when  expressed  in  the  x-y  system,  become 
i'  =  003(83  “  ")i  +  sin(83  ~  ir)j 

or 

i'  *  -cos  B3I  -  sinBjj  IB- 16) 

and 

=  sin  83I  "  cosBj  CB-17 ) 
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Fuze  Body  Configuration  No.  2 

Figure  B-4  is  used  to  define  the  angles  y^ ,  6^,  and  B^  of  configuration  no.  2. 
Angles  y 

Since  these  angles  are  defined  in  the  clockwise  direction  witli  respect  to 
the  body  fixed  X-  axis,  their  value  must  be  negative,  as  is  the  case  for  the  three 
pass  configuration  no.  2  of  appendix  A. 

According  to  equation  A-24  (app  A) 


Y 2  ~  -cos 


+  ^2  ~  (Rpl 
2  fli  fl2 


(B-iS) 


Angle  yj  is  given  by  equation  B-2, 


and  with  that 


y  becomes 


Y  3  =  Y  2  "  Y 3  CB-19) 

The  negative  sign  for  Y3  is  necessary  in  order  to  make  y3  negative. 

The  angle  y4  is  given  by  equation  B-4.  It  allows  the  following  expression  for 
angle  y4 


y4  =  y 3  -  y4  ( b-20) 

Again,  y£  must  have  a  negative  sign  to  make  y4  negative. 

Angles  6^ 

The  angles  62,  <53,  and  64  can  be  taken  directly  from  equations  B-6,  B-7, 
and  B-8,  respectively. 

Angles  8^ 

As  in  equation  A-28  (app  A) 


0t  =  7T  +  52 


CB-21 ) 


Further,  according  to  equation  A-29  (app  A) 


g2  =  y2  +  It  +  63 


( B-22 ) 
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Finally,  according  to  equation  A-30  (app  A) 


63  =  Y3  +  n  +  ^4  (B-23) 

Figure  B-4  shows  the  angle  y'  between  the  positive  x'-  axis  and  the  unit 
vector  n4.  This  angle  is  given  %y 


r 


p 


where 


it  +  e 


(B-24) 


e  =  n  -  64  - 


(B-25) 


Therefore , 


Y'  =  2n  -  S4  -  y4  (B-26) 

The  unit  vector  n4  can  again  be  expressed  in  terms  of  the  primed  coordinate  system 
as 


n4  =  cosYpi'  +  sinYpj  "  (B-27) 

In  p  manner  parallel  to  that  used  for  equations  A-3b  and  A-37  (app  A),  the  units 
vectors  i'  and  j  "  in  the  X-Y  system  are  expressed  as  fellows: 


i"  =  -cos  ^3!  -  sinBjj 


CB-28) 


and 


j'  =  sin83i  -  cos  B^j 


(B-29 ) 


Common  Computational  Expressions  for  Both  Configurations 

In  order  to  find  common  programming  expressions  for  all  angles  of  both  configu¬ 
rations,  with  the  exception  of  the  angles  y"  of  equations  B-14  and  B-26,  the  follow¬ 
ing  signum  function  s ^  is  again  introduced  p 


s6  -  +1  for  configuration  no. 


(B-30) 
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s  6 


-1  £nr  conf iguration  no,  2 


CB-31) 


This  lends  to  the  following  expressions: 


according  to  equations  B-l  and  B-18 


Y 2  M  s6cos 


-i 


a 


2 

R  2  “ 


(R 


J2l 


2  <Ri  R2 


JLLal 


(B-32) 


according  to  equation  B-2 


y  3 


cos 


«3-  (»p2  +  ro3)2 
2  fli  R  q 


(B-33) 


according  to  equations  B-3  and  B-19 


y  3  =  y  2  +  sgY3 


(B-34) 


according  to  equation  E-4 


Y 4  3  cos 


2  2  2 
R, +  «a-a 


ttt3  ri 


(B-35) 


according  to  equations  8-5  and  B-20 


Y4  "  Y3  +  S  gY 4 


(B-36) 


according  to  equations  B-6  through  B-8,  respectively 


62  *  CO 8 


(Rpl  +  rp2)2+ 

2^1  +  V 


(B-3 7) 


tS  3  =  cos 


&u  -  COS 


<Rp2  +  'p3>  2  +  *2-  a3 

2ir7(%rrf^> 

2 


2  iC2 2 

a  +  « , 


fl 


2a  3 


(B-38) 


(li-39 ) 


according  to  equations  B-9  and  B-21 


Pi  =  "  ~  s662 


(B-40) 


according  to  equations  B-10  and  B-22 


02  =  y2  +  it  -  s  663 


(B-41  ) 


according  to  equations  B-ll  and  B-23 


S3  =  y-z  +  *  -  sb<54 


(B-42) 


The  signum  function  s b  will  also  be  useful  in  conditional  statements  which  distin¬ 
guish  between  the  y"s  of  configuration  no.  1  and  2,  as  given  by  equations  B-14  and 
B-26,  respectively.^ 

The  expressions  for  the  unit  vectors  i "  and  j  "  are  identical  for  botli  configu¬ 
rations  according  to  equations  B-16  and  B—  1 7  as  well  as  equations  B-28  and  B-29. 
Thus , 


i'  *  -cos83i  ~  sinBjj 


( B-43 ) 


and 


j  '  »  sinB3i  -  cos  Bjj 


( B-44  ) 
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DYNAMICS  OF  PALLET  IN  COUPLED  MOTION 
(APPLICABLE  TO  BOTH  CONFIGURATIONS  WITH  APPROPRIATE  s6> 


Equation  A-95,  appendix  A,  which  gives  the  contact  force  Pn  between  pallet  and 
escape  wheel,  may  be  readily  adapted  to  the  present  two  pass  mechanism.  It  is  only 
required  that,  depending  on  configuration,  the  angle  y'  is  computed  according  to 
either  equation  B-14  or  equation  B-26.  Further,  the  distance  from  the  spin  axis  to 
the  pallet  pivot  is  now  given  by  <R4  rather  than  by  <R  5.  Thus  equation  A-9b  be¬ 
comes 


Pn  =  A^[1PRU^  +  (A21Z  +  IPRV)*  +  2  A2oZ*  +  “  A] 

“  mprcp  fl4<*>  sin(Yp-<|r^c)] 


lB-45) 


The  parameters  Aj  to  Ay,  AI8  to  A2l  and  I  may  essentially  be  taken  directly  from 
equations  A-70  through  A-77  and  equations  A-84  through  A-90,  appendix  A.  In  equa¬ 
tions  A-71  and  A-75  it  is  necessary  to  substitute  <R4  for  <R 


DYNAMICS  OF  ESCAPE  WHEEL  IN  COUPLED  MOTION 
(APPLICABLE  TO  BOTH  CONFIGURATIONS  WITH  THE  APPROPRIATE  s6) 

The  analysis  of  the  escape  wheel  can  be  taken  over  from  appendix  A,  once  the 
appropriate  signum  function  s^  has  been  chosen  for  the  given  configuration  and  the 
following  changes  in  nomenclature  have  been  made: 


*4 

becomes 

®  3 

becomes 

°2 

m4 

becomes 

ra3 

«4 

becomes 

«3 

becomes 

X3 

>4 

becomes 

^  3 

3  3 

becomes 

8  2 
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becomes 


becomes 


n  34  becomes 


n  23 


n.,_,  becomes 

N34 


N23 


F 34  becomes 


23 


x4 


becomes 


x3 


becomes 


P  3 


y4 


becomes 


y3 


b4 


becomes 


b3 


becomes 


Note  that  the  meanings  of  and  remain  unchanged  (eq  A-59  and  A-60). 

Similar  to  equation  A-117,  the  contact  force  equation,  as  derived  for  the 
escape  wheel,  then  has  the  following  form: 


-I  3<j)  +  F23A  15  "  T  3A  16 


P  = 
n 


1  7 


(.  H-46 ) 


where  now 


.  15  =  rfe3  -  M  [s  2(d  f  a  2)  +  P3<a  10  +  A  13>  J 


(B-47) 


A  16  *  M0j(A11  +  A  14) 


(B-48) 
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Aj  7  =  Aj  +  M  i  S4  3 1  +  pp  3(  Ag  +  A 1 2  ) 


and  similarly 


Ag  - 


-(Pi  s4-p  )sin(<f>-a-^3  )  +  ( l+pp 1  s4  )cos  (  4>—  oH- B 3) 

2 

1  +  P 


HO 


-p  (  1-  S2  )sin(02+®2  )  ~  (  1 +U  S2  ) cos (82+62  ) 


l  +  P 


AU  = 


siny  3  -  y  cosy  3 
2 

1  +  P 


A12  = 


(l+Pllis4  )sin(<J>-a-ifi3  )  +  (s4p  ^-p  )cos(<J>-a+£3  ) 


1  +  M 


13 


-  ( 1 +M  s2  )sin(B2+62 )  +  W ( 1~S2 )cos(B2+^2 ) 

2 

1  +  U 


U4 


-C0SY3  -  p slny 3 
2 

l  +  P 


(  B-49) 


(B-50) 


(B-51) 


(B-52) 


(B-53) 


(B-54) 


( B — 5  5 ) 
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COMBINED  COUPLED  MOTION  DIFFERENTIAL  EQUATION  FOR  ESCAPE  WHEEL  AND  PALLET 
(APPLICABLE  TO  BOTH  CONFIGURATIONS  WITH  APPROPRIATE  s6) 

Equations  B-45  and  B-46,  are  now  set  equal  to  each  other  in  order  to  obtain  a 
combined  differential  equation  for  the  escapement.  This  leads  to 


£  1  IS  +  LpRA  1  7U]  <f>  +  [XprA  .  7V  +  A  17A21Z  ]  $ 


+  2  A  17a2U7-^ 


=  a1Sa18f2J“  a  1  tA  1  8t  j  “  A  17w  r  A  ly  -  raprCp  qSin(  y'~  ijr-  (Pc)  "j 


(B-56) 


Again,  it  is  necessary  now  to  determine  an  expression  for  the  contact  force  p£ y 


DYNAMICS  OF  ROTOR  (Gear  No.  1) 

(APPLICABLE  TO  BOTH  CONFIGURATIONS  WITH  APPROPRIATE  s&) 


Figure  B-5  shows  a  free  body  diagram  of  the  rotor  of  configuration  no.  1,  which 
moves  in  a  counterclockwise  direction.  The  acceleration  of  its  center  of  mass  is 
given  by  equation  A-I22,  appendix  A. 

Since  the  motion  of  the  rotor  must  be  expressed  again  in  terms  of  the  escape 
wheel  variables,  use  must  be  made  of 


■f  l  =  1  $ 


(B-57) 


’1  ~  "31 


N31  * 


(B-58) 


where 


NP3  X  NP2 
3  1  *  n  x  N 

G2  G1 


(B-M) 
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The  rotor  angle  <^c  +  ^  is  expressed  in  terms  of  the  total  escape  wheel 
rotation  ^ 


1C 


+  =  *10  +  Njl  *T 


(B-60) 


(For  details,  see  program  description.) 


Newton's  force  equation  may  now  be  written  in  the  form  of  equation  A-128.  The 
tooth  forces 


F  1 1  =*  “F  i  &  !  2 


(B-61) 


and 


Ff2i  =  ns  ^2^2 


(B-62) 


have  the  appropriate  directions  to  account  for  friction. 


Thus 


-F  12n  12  +  US  jF  I2«n12  +  ^i1  "  bFyli  +  Fylj  +  pF^j 


2  _  2  _  _ 
mi  {  -  id  6^1  -  (w+  N31^)  rcl  [cos(  ♦lc+N31*r)±  +  sin(  <!>1C+N  31^)3  ] 


+  N 3  [-sln(<|>lr+N31<fer)i  +  cos C  4>j  r.+N 3 i^) J  J } 


*1 C  J 1  rT 


1C  3iTT' 


(B-63) 


where 


n  12  3  -sinCB!  -  6i)i  +  cos (3i  -  0 1) j 


nN12  =  _C08^1  "  8lH  "  sin(Bi  "  ® i>J 


(B-64) 

(B-65) 
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The  moment  equation  must  be  written  in  the  manner  of  equation  A-61  with  respect 
to  the  accelerated  point  0^.  The  pivot  reactions  F  ^  and  F  ^  are  written  such  that 
the  associated  friction  moments  retard  the  counterclockwise^  rotation  of  the  rotor. 


the  associated  friction 
This  leads  to 


~F  12*3,3  +  vejaiFia-  upi(fx1  +  Fyl> 


“>i<*>  flirclsin(<t,ic:+N3iV  +  1  iN 3 


(B-66) 


The  forces  F  and  F  are  obtained  after  the  simultaneous  solution  of  the 
following  component1  expressions  of  equation  B-63  for  F^  and  F^: 


Fl2sin(  Bj— 9 1)  +  gs  jF  12cos(  BrBi)  +  F^  -  uFyl 


*-  m 

■  m  -0)  -  (wFN31<())  rclcos(  <t>lc+N  3!^)  -  N31(j,rclsin(i)>lc+N3llj 


(B-67) 


-F  12cos(  63-63)  -  us  ^3 2810(63- 63)  +  Fyl  +  uF^^ 


m3  ^-(ort-Njil)  rclsin( ♦1C+’N 3 i<Jir)  +  N3i<jirC]lcos(tlc+N3i-{.T)J 


(B-68) 


Simultaneous  solution  of  the  above  furnishes 


2  ,  .2 

Fxl  “  F12A23  *  u  A24  ±  2  N  33IA25  ^  (N3ll)  A  25  ±  N  33^26 


( B-69 ) 


2  *■  2  ** 

F„i  "  F  12^27  ±  '■>  A28  ±  2  777  N31|a29±  (N33I)  A  23  ±  N  33^30 


(B-70) 


where 


Equations  B-69  and  B-70  are  now  substituted  into  the  moment  equation  B-66 


~Fl2Rbl  +  m  13  lF  12  "  MPl  fF  12^2^27)  *  «o  (A24+A26) 

2 

2  •• 

*  2  fuF  N3ll(A25+-A2<p  ±  (H31*)  (A254-A29)  ±  N31i)>(A26+A3o)  ] 


2  •• 
-  m  i<j>  M^sinU^+Na^)  +  1 3N 3 j d» 


(B-79) 


In  order  to  have  the  pivot  friction  moment  negative,  the  following  terms  on  the 
left  hand  side  of  equation  B-79  must  be  negative: 


UPlF  l2  (A23  +  A27) 


2 

UPlW  (A2i*  +  A28) 


MPl  (N3it)  (A25  +  A29) 


2 

2  1 1*1  P  1  n31  I  (A25  +  A2g) 


(With  N 3 a 4>  positive,  y  must  be  absolute.  See  discussion  following  equation  A- 
147).)  The  term  mPiNjj#  (A26  +  A30)  is  treated  in  the  same  manner  as  was  shown  in 
connection  with  equation  A-149,  appendix  A.  (See  the  discussion  in  appendix  F.) 
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The  above  considerations  give  the  moment  equation  the  following  form: 


“F  12  [Rbi  "  IS  ia  i  +  UP  l(A  2  J*"A  27>  ]  “  «*  UP  l(A24+A28) 

2 

2upjw  N  3  #  2 

- __  (A  25^A29)  ~  UPl  < N  3 1  <>)  (A  25+  A29) 


2  _ 

-  Him  ^r^sint^-H^i^)  +  (1 1±UPi(A2^A30))N31^  (B-80) 


Finally,  the  above  expression  is  solved  for  F12,  i*e.. 


..  2 


12 


"A 32  “A33N31$  ~  A3j,(N  31$)  “  A35sin(  <^C+N  3  j ij^)  -  IlRN31<f> 


31 


where 


a31 


A32 


A  33 


a34 


A  35 


Rbl  “  l«lal+  UPl(A23+  A  27) 
2 

UPi^  (A 24  +  A28) 

2 

2  |  ^  I p  1  jjjj’  (A  25  +  A29> 

UP  l(A  25  +  A  29) 

™1  fllrC1“2 


(B-81) 

CB-82) 

(B-83) 

(B-84) 

(B-85) 

(B-86) 
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(B-87) 


I 


1R 


1 1  +  |  u|pi(A26  +  A3q),  when 


♦ 


and  $  have  the  same  signs 


I1R  *  1 1  ~  |  u|pi(A26  +  A30),  when 


and  $  have  opposite  signs 


DYNAMICS  OF  GEAR  AND  PINION  SET  NO.  2 
(APPLICABLE  TO  BOTH  CONFIGURATIONS  WITH  APPROPRIATE  s6) 

The  free  body  diagram  of  gear  and  pinion  set  no.  2  of  configuration 
shown  in  figure  B-6.  Its  force  equation  is  written  as  follows  with  the  I 
force 


2 

T2  *  m2  (l2b3 


Thus 


*F23P23  ”  S2PF2?1N23  +  F12»  12  ~  »«  IF  +  Fx21  +  *jFy2i 


+  Fy2J  "  +  T2(cosy2I  +  sinyj)  -  0 


where 


n23  ■  sin(02  +  62)I  -  cos (  02  +  02)j 


n^j  -  cos(  82  +  62)I  +  sin(82  +  02)j 


(B-88) 


no.  1  is 
'Alembert 


(B-89) 


(B-90) 


(B-91) 


(B-92) 
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The  unit  vectors  nj2  and  n^2  were  6*ven  by  equations  B-64  and  B-65.  Hence 

F  32  “  “F23  n23  (B-93) 

Ff32  *  -32^23^23  (B-94) 


To  write  the  moment  equation,  let 

••  •• 

$2  “  N j2<J> 


where 


'32 


P3 


G2 


Then 

F 23^b2  -  VB  2F  23a  2-  F  12rb2  +  IB  lF  I2(d  ra  1)  +  W>2(Fx2+Fy2) 


(B-95) 


(B-96) 


"  l2N32<t> 


(B-97) 


The  bearing  forces  F^2  and  F  _  are  obtained  from  the  simultaneous  solution  of 
the  following  component  equations  ^or  the  force  expression  B-90: 


-F238in(  &2*Q2)  -  F23ys  2COs(  B2f02)  -  F  l2sin(  0  j-0  j) 


+  pa  iF  12cos(  B  i~ 0 1)  +  Fx2  +  jfy2  +  T  2COS y2  -  0 


(B-98) 


and 
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+  MS  12sin(  ^  1~ 6 +  Fy2  "  t*’x2  +  T23inY2  “  0 


(B-99) 


This  leads  to 


Fx2  *  F23A3b+  f12A37  +  T2a38 


(B-10U) 


Fy2  “  F 23  a39  +  f12  a40  +  T 2  A41 


(B-1U1) 


(all  signs  are  left  positive  in  order  to  furnish  positive  friction  moments  in  equa¬ 
tion  B-97)  where 


A 36  3 


A 37  _ 


A  38  = 


A  39 


l40 


A4  1 


(  1-p  s  2)sin(  B2+92^  +  h(  1+8  2)co8<  82+®2^ 


cos y 2  ~  iBinY2 


1  +  u 


p(l+s  2)sin(  92)  ~  (1“U  s  2)cos(  0zfe2) 


PCOSY2  +  sinY2 


1  +  M 


(B-102) 


(B-103) 


(B-104) 


(B-105) 


(B-lOb) 


(B-107) 


'  I 


101 


J 


tm 


Substitution  of  equations  B-100  and  B-101  into  the  moment  equation  B-97  and 
subsequent  solution  for  the  contact  force  F23  gives 


F  12  A42  “  T2a43  +  1  2  N32  <t> 


23 


l44 


(B-108) 


where 


A  1,2  =*  ~  U  [3  l(d  1  -  a  ])  +  P2(A37  +  A4o)  ]  (B-109) 

A43  =  up2(A3g+  A  43)  C  B  110) 

A44  =  Rfe2  “  U  ^8  2^2  ~  P  2^a  36  +  A  39)  j  (.B-lll) 


DYNAMICS  OF  COMBINED  SYSTEM  IN  COUPLED  MOTION 
(APPLICABLE  TO  BOTH  CONFIGURATIONS  WITH  APPROPRIATE  s6) 


In  order  to  obtain  a  single  differential  equation  for  the  total  system  in 
coupled  motion,  in  terms  of  the  escape  wheel  angle  <)>,  an  appropriate  expression  for 
the  contact  force  F23,  which  also  contains  the  contribution  of  the  rotor,  must  be 
substituted  into  equation  B-56. 

Thus,  first  substitute  equation  B-81  for  F12  into  equation  B-108,  which  is  the 
above  expression  for  F23: 


F  23 


-A  <♦  2A  3  2 

~*i r ' T2A'3 


A‘*Zfi3  3  .  A42A34  2  .2 

--T7TN31^  TIT N31* 


A42A35 

~A 7T 


1 42 


sin(^c+N31ir)  +  (I2N32  -  X77  I1rN3i)4* 


(B-112) 
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Now  equation  B— 112  is  substituted  into  equation  B-56,  and  the  final 
equation  of  coupled  motion  results 


••  2 

•  *-  • 

A45$  +  A46$  +  A47$ 


=  a48  -  Altq^in(((>lc+N31(hr)  +  A50pin(Yp~4r-<J>c) 


where 


a15^18  a42 

a45  -  ^  3A 1 8  +  ^1^  a^“  (I2N32  “  ^  riRN31> 


2  a15a18a34a42  2 

a46  =■  IprA17V  +  a17A2iU  +  *7^  N31 


A  ^3^42 

a-7  -  2  7^]-  aitA2oU  +  177^  n^i 


A15A18A32A42  A15A18A43  2 

A,*«  - - ^  T2  “  A16A:8T3  -  co  A17A19 


A49 


A  15A  laA  3^42 
A  3  1A  4  4 


2 

A  50  ”  mprCP  <R4A17“ 


differential 


(B-113) 


(B-114) 


(E-115) 


IB-llb) 


(B-117) 


(B-118) 


(B-119) 
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CONTACT  FORCE  EQUATIONS  FOR  COUPLED  MOTION 
(APPLICABLE  TO  BOTH  CONFIGURATIONS  WITH  APPROPRIATE  s6) 


The  contact  force  F23  is  given  by  equation  B-112 


23 


L  3  1 


_A32A42 


1 3 1 


“  T  2A 43  - 


A33^42 


A4 2 A  34 


- -  N  3 1 4>  -  — r - 

A  3 1  J  A  3 1 


2  .2 

N  3  1^ 


A  35A  42  A  42 

- IJT  sin(*lC+N3lV  +  (12N32  -  XJ7  I1RN31)* 


The  contact  force 

F  23  A 44  + 

F  12  - - 


is  found  with  the  help  of  equation  B-108 

T  2  A  43  "  F  2  N  32  4> 
a42 


The  contact  force  Pn  is  given  by  equation  B-46 


P 

n 


-I3  $  +  F 23  A15  -  T3  A16 


1 1  7 


If  this  force  is  desired  in  terms  of  the  pallet  variable  i|>,  equation  A 
used  with  the  understanding  that  it  now  has  the  following  form: 


PR 


.2  m  2 

♦  +  A21ip  +  2  A 20^  +  0)  A  19 


-  mprcp  u)  sin(  yp-iM'c) 


A  18 


where  <R^  is  used,  together  with  y'  as  defined  earlier  in  this  appendix, 
previously,  all  other  parameters  may  be  taken  from  appendix  A. 


(B-120) 


(B-121) 


(B-122) 


91  must  be 


( B-123  ) 


As  stated 
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DIFFERENTIAL  EQUATIONS  FOR  FREE  MOTION 
(APPLICABLE  TO  BOTH  CONFIGURATIONS  WITH  APPROPRIATE  sfe) 


Free  Motion  of  Pallet 

Let  P  =  0  in  equation  B-123.  After  some  rearrangement,  the  free  motion 

differential  equation  for  the  pallet  becomes 


.2 

A  5 1 4*  +  A2i<|>  +  A52<|/  =  -ASj  +  As^sinty'-iJ'-'l^) 


(B-124) 


where 


A51  =  1 PR 


(b-1251 


=  2  lu)|  A20 


(B-12b) 


,\5J  -  a)  A  1;) 


(B-127; 


A  54  =  “  rap  rCF  6'4 


(B-12B) 


(Again,  the  angle  Yp  is  defined  for  both  configurations  by  eq  B-14  and  B-2b.) 


Free  Motion  of  the  Escape  Wheel,  Gear  Train,  Rotor  System 

By  setting  Pn  ■  0  in  equation  B-46,  the  differential  equation  of  the  system 
without  the  pallet  is  obtained.  This  leads  to 


x3  ♦  *  A 1 5  E  23  *  T3  Aie 


(B-l 29) 


Equation  B— 112  is  then  substituted  into  the  above  for  F23»  The  resulting 
differential  equation  is  given  by 


(B-130) 


A55 <t>  +  As6|  +  A57J  ■  Asfi  “  A598in(^c  +  N 3 j 


where 


A  1 5  A15\42 

A55  "  X3  "  N32l2+  X7]a^  N31I1R 

A  1 5A  34A  42  2 
A56  '  ~ *7^-  "31 


A  1&A  33A42 

A57  '  ”A3'l‘~  "3‘ 


A15A43  A  15A  32A42 
A58  *  -auT3  -  Tj  -  ■  *7,*-— 

A15A35A42 

■  "x3>-^r 


(B-131) 


(B-132) 


(B-133) 


(B-134) 


(B-135) 


CONTACT  FORCE  EQUATIONS  FOR  FREE  MOTION 
(APPLICABLE  TO  BOTH  CONFIGURATIONS  WITH  APPROPRIATE  s6) 


known 


Equation  B-129  may  be  solved  for  the  free  motion  contact  force  Fp23*  once  $  is 


I3  <j>+T3Al6 


(B-136) 


Equation  B-108,  which  was  derived  from  the  force  and  moment  equations  of  gear 
and  pinion  set  no.  2,  may  be  modified  to  obtain  the  free  motion  contact  force  Fp^ 


A44  +  T  2  A43  -  I  2  N32 


(B-137) 
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It  must  be  understood  that,  both  In  equations  B-136  and  B-137,  the  angular  acceler¬ 
ation  of  free  motion  <j>,  as  obtained  from  the  solution  of  equation  B-130,  Is  used. 


CHANGES  IN  IMPACT  EXPRESSIONS 
(APPLICABLE  TO  BOTH  CONFIGURATIONS) 


Both  the  rotor  and  gear  and  pinion  set  no.  2  must  be  referred  to  the  escape 
wheel  shaft.  Thus,  similar  to  equation  A-212,  appendix  A 

iST0T  “  I3+  *2  «32+  N31  (B-138) 


See  equations  B-59  and  B-96  for  the  above  transmission  ratios. 
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DIMENSION  PHI (2 ) ,  DPMI ( 2 ) ,  PftMT<5) 
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1BOA,DELrA,PHITOT,PHIPR,N4l , N42.N43, OMEGA .OM2.RC1 , PH1 1 C , T EST 1 , T EST2 
2 , TEST  3 ,NG1 , NG2,NG3,NP2,NP3,NP4,CAPR31  .  Q\PRR2 ,CaPRB3 , RB2 , RB3, RB4 , TH 
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COMPUTATION  OF  A1 .  A2,  AND  A3 


Computer  program  SANDA3  (cont) 
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Computer  program  SANDA3  (cont) 
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Computer  program  SANDA3  (cont) 
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Computer  program  SANDA3  (cont) 
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PrtECKDlNG  PACK  BUNK -NOT  FILMED 


Since  nil  friction  forces  of  the  escape  wheel  to  rotor  gear  train  reverse 
direction  during  a  possible  motion  reversal  following  an  impact,  the  introduction 
of  a  negative  coefficient  of  friction,  i.e.,  a  -p,  will  satisfy  for  computational 
purposes  the  associated  force  and  moment  equilibrium  expressions.  (This  method 
need  not  be  resorted  to  for  the  pallet-escape  wheel  interface,  since  there  motion 
reversal  is  accounted  for  by  the  signum  function  s^.) 

The  above  will  be  proven  to  hold  for  both  the  tooth  contact  and  the  pivot 
friction. 


Tooth  Contact  Friction 

The  direction  of  the  friction  force  of  the  pinion  on  a  gear  of  a  mesh  is  the 
same  as  that  of  the  relative  velocity  of  the  contact  point  on  the  pinion  with 
respect  to  the  contact  point  on  the  gear  (ref  E-l).  This  relative  velocity  is 
obtained  by  way  of  the  vector  equation 


C  /C 
P  g 


+  V, 


(E-l) 


where  (fig.  E-l) 


C  /C 
P  g 


absolute  velocity  of  contact  point  Cp  of  pinion  tooth,  with 
direction  normal  to  line  connecting  pivot  point  0  and  contact 
point  C  (contact  points  C  and  CT  are  coincident; 

r  r  6 

absolute  velocity  of  contact  point  C  of  gear  tooth  with 
direction  normal  to  the  line  connecting  the  pivot  point  Og 
and  the  contact  point  C_ 

O 

relative  velocity  between  contact  points  Cp  and  Cg.  The 
direction  of  this  velocity  is  normal  to  the  line  of  action 


Velocity  triangles  according  to  equation  E-l  for  contact  of  a  simple  mesh 
above  and  below  the  pitch  point  for  two  cases  is  shown  in  figure  E-l.  In  figure 
E-l (a),  the  gear  is  the  driver.  In  figure  E-l(b),  where  the  motion  is  reversed, 
the  pinion  becomes  the  driver.  Figures  E-l (a)  and  E-l(b)  differ  with  respect  to 
the  direction  of  the  relative  motion,  and  with  that,  the  direction  of  the  contact 
friction  force  for  comparable  contact.  This  difference  in  the  direction  of  the 
friction  force  may  be  expressed  by  reversing  the  sign  of  the  friction  force.  For 
computer  purposes,  this  is  best  accomplished  by  making  p  positive  when  the  gear 
is  the  driver  and  making  it  negative  when  the  pinion  becomes  the  driver.  The 
normal  force  between  the  gear  and  the  pinion,  which  has  the  direction  of  the  line 
of  action,  is  not  dependent  on  the  driver. 
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Direction  of  rotation 
of  driver 


Pivot  Friction 


The  free  body  diagram  of  the  combined  gear  and  pinion  set  no.  2  in  clockwise 
rotation  is  shown  in  figure  E-2.  The  associated  equilibrium  equations  B-39  and 
B-60,  which  are  used  here  as  a  general  example,  may  be  rewritten  in  the  following 
manner : 


Fx  +  =  A1  + 


pF  12COSI 


and 


(E-2) 


-uFx  +  Fy  =  A2  +  uF23SIN2  -  pF  US1N1  (E-3) 

where  Aj  and  A2  stand  for  the  normal  tooth  contact  and  pivot  forces  which  do  not 
contain  the  coefficients  of  friction  and  do  not  change  direction  on  motion  rever¬ 
sal.*  The  trigonometric  functions  refer  to  the  associated  angles  gj  -  0 1  and 
82  +  02*  The  signum  functions  s^  and  s2  are  omitted  for  the  present  discussion 
since  their  influence  is  not  relevant  in  this  context.  They  change  the  signs  of 
all  friction  terms  on  the  right  sides  of  equations  E-2  and  E-3,  depending  on 
whether  contact  is  made  before  or  after  the  pitch  point. 

Assume  now  that  the  direction  of  the  gear  set  is  changed  to  be  counterclock¬ 
wise.  This  causes  the  tooth  contact  friction  forces  (as  shown  earlier)  to  be 
reversed  in  sign.  Furthermore,  the  two  pivot  friction  forces  are  also  reversed 
in  sign.  The  above  force  equilibrium  equations  now  have  changed  signs  in  all 
friction  terms,  i.e., 


pF  +  F  =  A2  -  pF  23SIN2  +  pF12SItn  (E-3) 

X  y 

Clearly,  this  effect  of  the  change  of  the  direction  of  all  friction  forces  can  be 
expressed  by  making  the  coefficient  of  friction  negative  for  computational  pur¬ 
poses  . 


Equation  A-2,  reference  E-l  shows  that  if  the  assumptions  of  the 
directions  of  the  normal  pivot  forces  prove  wrong,  the  associated 
friction  forces  still  furnish  friction  moments  which  have  the 
appropriate  direction  of  rotation. 
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This  idea  also  holds  true  in  general  for  the  moment  equation  and  is  illus¬ 
trated  with  the  help  of  the  moment  equation  B-58  for  gear  and  pinion  set  no.  2, 
which  for  clockwise  rotation  is  of  the  following  general  form: 


B  -  pF23a2  +  pF  12(d  pa  +  up2(Fx2  + 


Fy2>  =  1&32* 


(E-b) 


where  B  is  similar  to  A^  and  A2  in  equations  E-2  and  E-3  by  representing  moment 
terms  which  do  not  contain  the  coefficient  of  friction,  and  whose  sign  does  not 
change  on  motion  reversal.  In  the  same  vein,  the  signum  functions  s^  and  s.,  are 
omitted . 

If  the  direction  of  rotation  is  counterclockwise,  the  above  expression  be¬ 
comes 


B  +  pF  2  3a  2  ~  Pf  12^  l_a  l) 


PP  2(f 


x2 


Fy2)  “  I  2*<324 


(E-7) 


Thus,  the  change  in  the  direction  of  the  friction  moments  can  also  be  ex¬ 
pressed  by  making  the  coefficient  of  friction  negative  for  computational  pur¬ 
poses. 


Finally,  the  results  of  these  new  force  and  moment  equations  for  the  various 
components  may  be  obtained  by  changing  the  signs  of  all  coefficients  of  friction 
in  the  results  of  the  original  expressions  which  were  derived  for  forward  motion 
of  the  train. 
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APPENDIX  F 

A  COMPUTATIONAL  RULE  FOR  THE  DETERMINATION  OF  THE  SIGN 
IN  THE  EFFECTIVE  MOMENT  OF  INERTIA  I 

iR 
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Whenever  the  center  of  mass  of  a  fuze  train  component,  such  as  a  pallet  or  a 
rotor,  does  not  coincide  with  its  pivot  axis,  there  will  be  a  pivot  friction 
moment  which  is  proportional  to  the  angular  acceleration  of  the  component.  Just 
like  any  other  friction  moment,  its  direction  must  be  opposite  to  the  angular 
velocity  of  such  a  rotating  part.  The  following  shows  how  the  appropriate  choice 
of  the  sign  in  the  expression  for  the  effective  moment  of  inertia  l^g*  of  the 
component  assures  that  this  friction  moment  will  always  retard  motion.  While 
this  discussion  holds  in  general,  it  makes  specific  use  of  the  pallet  moment 
equation  A-7H,  appendix  A.  This  expression  becomes  in  condensed  form 


M  ±  M  =  1  ip  C  F- 1 ) 

S  AA  p 

where 


Mg  =  the  sum  of  all  moments  on  the  pallet,  with  the  exception  of  the 
friction  moment  which  is  proportional  to  the  angular 
acceleration  of  the  pallet.  The  moment  Mg  is  sufficiently 
larger  than  the  moment  so  that  it  will  at  all  times  decide 

the  direction  of  the  angular  acceleration  >ji 


MAA  =  Pp^+ABH  =  A22'l' 


(F-2) 


i.e.,  A>22  is  always  positive.  Note  also  that  the  siguum  function 
s^  has  been  omitted. 


Figure  F-l,  a  through  d  represents  free  body  diagrams  of  a  simplified  pal¬ 
let,  or  rotor,  with  various  directions  of  M  ,  and  M  . 

j  /vA 

In  figure  F-la  the  moment  Mg  is  assumed  to  be  CCW,  i.e.,  positive.  The 
angular  velocity  i}i  is  also  assumed  to  be  positive  and  therefore  the  friction 
moment  must  be  negative.  Equation  F-l  becomes 


Ms  -  maa 


=  I 


(F-l) 


The  letter  i  In  the  subscript  stands  for  the  component  number  or  letter. 
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tHECSDUO  MCS  MJiMC-NOT  TJiJCD 


Ill  general,  the  sign  of  ip  is  never  predetermined  in  a  differential  equation  of 
this  type,  but  it  is  decided  by  the  sign  of  the  resultant  moment.  Ln  this  spe¬ 
cific  work  it  must  be  remembered  that  the  direction  of  ip  is  decided  by  Mg  alone 
and  that  the  friction  moment  must  act  opposite  to  the  direction  of  the  angu¬ 

lar  velocity 

is  substituted  into  the  above  expression  witli  a 


(K-0 


F-4  now  be  used  to  find  equivalent  expressions  for  the 
Thus,  from  equation  F-3 

U-'-b) 


(F-h) 


In  the  first  form  the  absolute  value  of  the  angular  acceleration  is  reduced,  from 
the  value  it  would  attain  due  to  the  moment  Mg  alone,  by  the  presence  o  f  Maa .  I n 
the  second  torm,  which  must  be  used  in  the  computations,  i);  is  reduced  i.i  absolute 
magnitude  by  the  equivalent  presence  of  a  larger  effective  moment  of  inertia  lpp 
=  Ip  +  (Tile  subscript  PR  is  used  for  the  pallet.) 

Now  consider  the  condition  of  figure  F-lb  where  Mg  is  also  positive, 
but  is  negative,  and  with  that  M^  must  be  positive.  This  leads  to  the  follow¬ 
ing  moment  equation: 

M„  +  M . .  =  I  i|»  U'W) 

S  AA  p 

or,  with  equation  F-2  and  a  positive  <|» 


If  equation  F-2 
positive  ip,  it  becomes 


Mg  =  <Klp  +  A22) 


Let  equations  F-3  and 
angular  acceleration  ip. 


♦ 


MS  "  haa 


and  from  equation  F-4 
Ms 

^  ~  I  +  A  22 

p 


M 


S 


(F-d) 


The  above  expressions  are  again  used  to  find  two  equivalent  forms  tor 
From  equation  F-7  one  obtains 


*  = 


Mc  +  a.. 
S  AA 


I 

P 


(F-9) 
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while  equation  F-8  gives 


■h  = 


(F-lu) 


tn  the  first  form  the  absolute  value  of  \j)  is  increased  by  the  added  presence  of 
M\\  over  that  due  to  Me;  alone.  The  second  form  shows  this  increase  in  ijj  by  way 


of  the  reduced  effective  moment  of  inertia  I 


PR 


"  1p  "  A22* 


In  figures  F-lc  and  F-ld  Mg  is  shown  to  act  in  a  clockwise,  i.e.,  negative 
direction,  and  with  that  iji  is  also  negative.  Consider  the  condition  of  figure  F- 
lc  first.  The  moment  equation  must  take  the  negative  into  account  by  way  of  a 
positive  M...  Thus,  from  the  free  body  diagram 


-  Mg  + 


AA 


=  V 


(  F-l  1  ) 


To  keep  all  signs  correct  on  substitution  of  equation  F-2  and  the  negative  \'i,  the 
above  expression  must  take  the  form 


-  Ms  ~  A  22 


1F-12) 


or 


-  Mg  =  -iK!  + 


IF- 13) 


The  absolute  value  of  ij>  is  then  given  by 


v  = 


3 

i  + 

p 


'22 


U'-lh) 


Again,  this  value  is  reduced  in  absolute  magnitude  by  the  increase  of  the  ettee- 
1 1 ve  moment  of  inertia  to  l  (J  ^  =  Ip  +  A^*  This  is  the  equivalent  to  a  decrease 
in  the  absolute  value  of  the  resultant  moment. 


Figure  F-ld  shows  the  pivot  friction  moment  in  the  same  direction  as  Mg 

since  ip  is  positive.  The  tnoim  nt  equation  becomes 


-Mg  -  M 


AA 


I  'P 


(MO 


l/pon  subs l  i  t ii t  i on  of 


equation  F-2  together  with  the  negative  \)j,  one  liuds 


1<J? 


i 


-Mc  +  A 22  (“<lO  =  I  (“♦) 

S  p 

or 

-Ms  =  -  <P(lp~A22) 


The  absolute  value  of  ip  now  becomes 


1 F- 1 o ) 


(F-1 7 ; 


<p  = 


(F-18) 


It  is  increased  by  the  decrease  in  the  effective  moment  of  inertia  to  lj,  ,  = 

Xp  ~  A22* 

These  four  cases  allow  the  following  computational  rule  for  the  determina¬ 
tion  of  the  sign  of  the  parameter  A22  in  the  effective  moment  of  inertia. 

From  the  equations  F-6  and  F-14,  together  with  the  associated  figures  i-la 
and  F-lc,  respectively 


PR 


=  I  +  A 


22> 


when  Ji  and  ip  have  identical  signs 


Further,  from  equations  F-10  and  F-18,  together  with  the  associated 
and  F-ld,  respectively 


tF-19) 

figures  l-'-lb 


I  =  I  -  A22,  when  tp  and  ip  have  opposite  signs 
PR  p 


(F-2U) 
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